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II. Operation : 

A) Load hex tapes A and B (See storage allocation A and B for tape 
contents). Note exception: 10. I4 is not on tape A or B. 

B) Insert sample problem data tape in the typewriter and depress 
break points ii and 8. Start the tape reading via a program 
input routine (10.0, 10.3 or 10. U). All start fill and trans- 
fers are on the tape. 

C) The computer will halt on the stop and transfer instruction 
(.0003000*). A start compute signal will initiate the com- 
putation (or B.P. 32 could have been depressed). 

D) The program will proceed to the completion of the regression 
analysis and then stop in location 1;006 (XZ0000). 

E) A start compute signal will cause an entry into the subroutine 
to compute the residuals, (Y-Y ca -^). Although the regression 
equation coefficients are stored correctly by the subroutine 
that computes them, they are in floating point and must be 
entered in the correct order in fixed point before the first 
record of the original data is read. This information is on 
the sample problem data tape. 



III. Format : 

A) Set carriage return stop at J4. 

B) Set tabs at 12, 20, 28, 36, ijlt, etc. n+1 tabs. (8 numbers apart), 

C) For very large problems the automatic C.R. must be used or all 
spacing tabs in program 1 and 3 removed (see write ups) and 

a column type printout used. 

D) Since break points separate every major phase of this group 
of subroutines, these spaces may be used for exits to heading 
printing operations or other calculations, without disrupting 
actual program steps. 

IV. General : 

A) The subroutine designated as "No. 2" (scaling) has not been used 
in this sample problem since then, in general, special handling 
would be required in the residual computation. 

B) The subroutines were assembled except for locations 3219-3223 
and 1|006 by filling the following subroutine Lo on top of the 
previous subroutine L|.. 
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FIXED POINT DATA 
COMPUMATIX, INC. 
SAMPLE PROBLEM STATISTICS SUBROUTINES 



j 0006232' 

xzOOlO' 

xz0003 ' 

xz0005 ' 

xz6200 ' 

xzl*200' 

xzl200' 

xzl+201' 

xziilOO 1 

xzOOOO 1 

xz6l32' 

xz2800' 

xzlO-32 1 



q of printout 

no. of variables 

no. of records 

Lo. of record 

Lo of A data matrix 

Lo. of floating point 

Lo. of means 

Lo. of standard deviations 

Lo. of sq. Rij matrix 

Lo. of scale factors 

Lo. of matrix inversion (29.0) 

Lo. of coefficients of regression equation 



depress break points h and 8 for straight thru operation 
to the residual computation 

.0003000' 
data 



0+0^6200 ' 1 " k * 1 ' -0000000 ' ' 
0+0^6200 ' 1 ' 1 ' k ' 3 ' -0000000 ' ' 
0+056200 ' 1 ' 3 • 3 ' 2 ' -0000000 ■ ' 
0+056200 ' 1 ' 6 ' 2 ' 5 ' -0000000 ' » 
0+056200 ' 1 ' 8 " h ' -0000000 ' ' 
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N 

5 . 00000 



ZXi 

18.00000 
? 

AA-, 

110.00000 



zx 2 

13.00000 
25.00000 

U5. ooooo 





ZY 


15. 


,00000 




ZX ] _I 


71. 


.00000 




ZX 2 Y 


32, 


.00000 



ZY^ 

, 55.ooooo 



x l 

.3600000 


01 


x 2 
. 2600000 


01 


Y 
.3000000 


01 




<TXy 

.3006659 


01 


.11*96663 


01 


6 1 
.11*11*213 


01 




r ll 
.1000000 


01 


r 12 

.9688981- 


00 


r iy 
.7996129 


00 








r 22 
.1000000 


01 


r 2y 
.66111376- 

r 

yy 
.1000000 


00 
01 








Inver 


se of 


r • - matrix 






.6081*1*1*3 


02 


.1*7591*27 


02 


.1717131*- 


02 




.1*7591*27 


02 


.3900739 


02 


.1225603- 


02 




.1717131*- 


02 


.1225603- 


02 


.6623821* 


01 




Pxl 




b xl 




r iy,2 




^bxl 


.2592361 


01 


.121931*1* 


01 


.85531*12 


00 


.5222230 


Px2 




b x2 




r 2y,2 




<5"bx2 


.1850295 


01 


.171+836U 

b 

.5935385- 


01 
01 


. 7621*686 


00 


.101*9098 

<*~b0 

.38851*88 


R 
.921U280 


00 


0~EST 
.51*91*909 


00 


R 
.8351*996 


00 


A 

<5~EST 
.8688215 



00 



01 



00 



00 
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regression coefficients 

6+051+132 ' -5935385 ' 12193UU * 17U8365 ' -0000000 ' 
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data repeated to compute residuals 
printout order Y, Y calj (Y-Y cal ) 

0+056200 ' 1 ' ' k ' 1 ' -0000000 ■ ' 
0+056200 ' 1 • 1 ' k ' 3 ' -0000000 ' ' 
0+056200 , l , 3'3 , 2 , -0000000" 

0+056200 ' 1 ' 6 • 2 ' 5 ' -0000000 « ■ 
0+056200 • 1 ' 8 ' ■ k ' -0000000 ' • 



Y 


Y 
cal 


(Y-Y .,) 
v cal 


1.00000 


1.05808 


~ .05808 


3.00000 


2.277U2 


.72258 


2.00000 


2.9677U 


- .9677U 


5.00000 


U.877U1 


.12259 


h. 00000 


3.81937 


.1806U 
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3000-3218 



Fixed Point 
A Matrix Prep 
Print & Float 



3219-3263 

Matrix 
Scaling 



JC 



3300 



3515 



Calc. & Print 
\f^ r ij 
Floating Point 



3516 



3606 



Convert A r^4 
matrix to 

□ r. . matrix 

(Modified) 



3607 



3660 



Invert & Print 
inverse D r^ 
Matrix 



3662 



U005 



Calc. & Print 

Pxi > ^xi > r xiy 
Sbxi,bQ, R, S, 



r# 



I4007 1 ■ U060 



Calc. & Print 
(Fixed Point) 

l> Y C ai (y-y C ai) 



T 



End 



Statistics Subroutines 
Fixed Point Data 
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TITLE: 



AUTHOR : 



DATE: 



Least Squares Fixed Point Triangular Matrix Preparation 

William F. Burggrabe 
Compumatix, Incorporated 

August 31, 1961 



PURPOSE: Given N records of n variables of the form 1, X-^, Xg, ... X n , 
~ this subroutine will prepare a fixed point 

triangular matrix. The matrix rows are stored sequentially 
in the following format: 



Z X 
Z X 



Z X, 



Z X.]X 2 



Z X< 



Z X^ 



zx x x n 



Z X 



n 



The matrix is then put into floating point form by sub 25. OR 
following an optional (T.C. down - no print out) fixed point 
matrix print out. 



RESTRICTIONS: 



1) All data are assumed to be at the same q, summations are 
computed at twice the data q. 

2) The following subroutines are required: 

Linkage Location Subroutine Fixed Storage Location 
Lo + 155 and 156 11.2 0300-0563 

Lo + 56 and 57 12.1 (a) 0600-0832 (0850) 

Lo + 207 and 208 25. OR 0900-1163 

3) The following information must be supplied: 

Location : 

6232 2 x q of data 

6233 number of variables 
623U number of records 

6235 Lo record 

6236 Lo A matrix 
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1 

tl) If any summation exceeds the maximum possible number size 
for 2.q, an overflow stop will occur in Lo + 105, (in gen- 
eral the sum of squares terms will create the most diffi- 
culty) and the problem then must be rerun with a greater 
data q. 

5) No limit to the number of records provided no overflow is 
encountered; the number of variables is limited only by 
available machine storage. 

Required storage: Data Lo to Lo + n + 1 

»« -.- ■ t 4. t (n+l) (n+2) 
Matrix Lo to Lo + -i '- — i '- 



6) Data Format: 
Initial record. 

P+ qqRRRR ' 1 ' X-^ X 2 ' ••• V -0000000" 

Following records: 

P+ qqRRRR+1' 1^ X 2 ' ... X ' -0000000" 

where RRRR = Lo record as stored in (6235) 
qq + Half of 2q as stored in (6232) 

7) Only the record being processed is stored. 

GENERAL INFORMATION ON METHOD : 

Since little mathematical description need be given, a 
description of the program sequence will be outlined here: 

1) Initialization 

Clear the matrix storage area to zero 

2) Semi-optimum matrix preparation. 
N-records, n variables 

3) Optional print out of fixed point matrix 
(Transfer control up) 

k) Conversion of fixed point matrix to floating point. 
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A break point stop separates each phase of the operation, as 
well as the processing of each record. 



CODING INFORMATION: 



A) Storage: Subroutine 2 tracks 20 sectors 

External storage 6232 thru 6236 as under restriction 3 
and the following . . . 

6260 ctr 

6261 n + 1 @ 29 

6262 n + 2 @ 29 

6263 (n+1) (n+2) @29 

2 

Data Lo to Lo + n + 1 

Matrix Lo to Lo + (n+l) (n+2) 

2 

B) Linkage and calling sequence: Since this is one of a 
group of statistical subroutines, location 6232 through 
6263 have been reserved for initialization information 
and no calling sequence is required. 

Linkage R Lo + 219 

U Lo 

C) Input: As described under restriction 3 and 6 

D) Output: The computed matrix is stored sequentially- 
beginning in the specified location. It may be printed 
out in fixed point in the following format: 







as per 12.1a 




N 


z x x 


Z* Ap • • • 

as per 12.1a 


ZX n 


Tab 
Tab 


z xf 

Tab 


Zj An Ap • • • 

Zj Ap • • • 


zx 2 x n 



etc. 
(See sample problem) 



The tabs may be suppressed by changing Lo + 001+9 to 
U (Lo + 55). 
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E) Location of constants: 
Lo + 110 1 @ 3 



Lo + 11U 


1 @ 29 


Lo + 115 


1 @ 29 


Lo + 116 


3wwj 


Lo + 117 


3wwj 


Lo + 123 


1 @ 29 



F) The actual matrix preparation time is approximately: 

[ ( n+1 ) ] n S ec: compute time, + ( n ) N sec. 
12 3 

Data read, time (Photo reader). 

G) Program Stops: 

Lo + 3k (Break Point k) after initialization 

Lo + 37 (Break Point 8) after completion of matrix prep. 

Lo + l%k (Break Point k) after each record is processed. 

Lo + 200 (Break Point k) after print out of matrix 

Lo + 218 (Break Point 8) after floating matrix. 
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FIXED POINT MATRIX PREPARATION 



Operating Procedure: 



Prepare data tape as per instructions under input and fill the required 
storage locations in 6232-6236 as under restriction 3. Insert data tape 
in reader and transfer to routine Lo or use a main program with the re- 
quired linkage. 

After inialization a stop (BP h) in Lo + 3U will occur. Depression of 
start compute will cause the reading of the first record, the necessary- 
computations for that record and an exit to (BP k) in Lo + l5U prior to 
reading the next record. (This feature allows for machine stoppage on 
very large problems). 

After all records have been processed, the program will halt on BP 8 
in Lo + 37. If it is desired to eliminate the fixed point matrix print 
out, depress transfer control before continuing the program. This 
causes a transfer to Lo + 200 (BP k) the end of the print sequence. 
(NOTE: Depression of the transfer control switch during the printing 
operation will cause an exit from the print routine at the end of the 
particular line being printed ). 

A start compute signal will then cause the fixed point numbers to be 
floated. On completion of the float sequence, a terminal stop on BP 8 
in Lo + 218 is executed, and a start compute signal will transfer 
control back to the main program. 
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SAMPLE PROBLEM FIXED POINT A MATRIX PREP. 



POOL Program no. - F2-129 

1 



0006232'xzOOlO' 
0006233'xz0003' 
000623ii , xz0005 I 
0006235'xz6200' 
0006236 'xz^OOO' 



print q 
no. variables 
no. records 
Lo record 
Lo A matrix 



> External Storage 



.0003000' transfer to matrix prep, data follows 



0+056200 ' 1 " k ' 1 ' -0000000 • • 
0+056200 ' 1 < 1 ' k • 3 « -0000000 < ' 
0+056200 > 1 ' 3 ' 3 1 2 • -0000000 ' ' 
0+056200 • 1 ' 6 « 2 ' 5 ' -0000000 ' ' 
0+056200 ' 1 ' 8 " k ' -0000000 ' ' 



-\ 



Data 



J 



5.00000 



18.00000 
110.00000 



13.00000 
25.00000 
li5.ooooo 



15.00000 
71.00000 
32.00000 
55.00000 



TITLE: Fixed Point A Matrix 
Prep, and Float 



003i* 



0037 

8 






No 



0200 



P 



Initialize 
Sub . 
0000-0025 



Program No. - F2-129 -1 
August 31, 1959 



Clear Matrix 
Storage area 

0026-0033 



Read Record 
0155-0156 




Compute and 
store suns of 
cross products 
Reset 




0151* 



3P 



Matrix 
Prep. 



Down 



Print Matrix 
Row 

00Li9-0063 




0216 



BP 

# 8 



EXIT 
A 



Float Matrix 
0201-0217 



Yes 
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JOB NO. 


PROGRAM NO. 

F2-129 


PROGRAM PREPARED BY: 

Burggrabe 


PROGRAM CHECKED BY: 

POOL Review 


DATE M 

U-18-59 


PROBLEM: 

Given! N records of n variables; prepare A Fixed point matrix & Float. 


TRACK 

00 


PROGRAM INPUT CODES 




h 


LOCATION 


INSTRUCTION 


a. 


h 
w 


CONTENTS 
OF ADDRESS 


NOTES 


OPERATION ADDRESS 


1 


1 
















1 


X 








,0 l0 |0 


i i X|B | 6, 2,3 3 




n@ 29 


n~ 

(No. of variables) 






1 |0 i 1 


i i ,A 0, 1,1 ,1* 




1® 29 








1 1 ,2 


i , X,H 6, 2,6 p. 




(n+1) 29 








1 !<> I 3 


i i i A - 0,1,1 fa 




X 


1® 29 








1 |0 | 4 


, , X,H , 6, 2,6 ,2 




(n+2) @ 29 








1 1° 1 5 


i i X, N | 6 2 , 6 .1 




(n+1)® 29 


give 3 
(n+l)(n+2) @ 27 






1 l° l 6 


i i i M | °i 1 1 1 i° 




IP 3 e iv 


e(n+l)(n+2) a 0Q 






1 l° 1 "> 


i i x i c , 6| 2 ,6 i3 




X 


(n+l)(n+2) 


i 29 


n " 1 


1 1 I | 1 1 1 




1 1° I 8 


i i X, B . 6 2 , 3 5 




Lo Record 








1 1° 1 9 


1 1 1 Y 1 °l 1 1 ° 1° 






] 






1 l' |° 


. ( | T | ( 1 ( 1 






Mlnit. Matrix Prep.) J 






i i t , t 


i i i T - q l , l ,3 




X 




, 






1 1 1 | 2 


i | X, A. 6, 2,6,1 




(n+1) @ 29 - 


■*■ Record L„ + 1 






1 1 1 i 3 


I i i Y | Q 1,1 ,1 












i i 1 i 4 


i i X,B, 6 2,3 ,6 




Lo A Matrix 








1 M i 5 


i i i Yj Q 1,0 i5 




X 










l M l 6 


i i i T. q 1,0,6 






/(Init. Matrix Prep.) ' 






1 |1 i 7 


i 1 1 Y l q 1 , 1 , 8 






1 






1 |' i 8 


i i J, q 0,2,7 






(Init. Clear Drum) 






1 J 1 9 


i i i x , q 0,5,5 




X 




(Init. Print) 






1 |2 ,0 


i i , T. q 2, 0,6 






X Init. Float 






i |2 i 1 


i i i Y i Q 2 i 1 . 






J 






1 i 2 i 2 


i , X, A. 6, 2, 6,3 




(n+1) (n+2) _ 
*r — ' <3> 


29 






I 1 2 i 3 


! 1 \ C | 4 2 i 6 ,° 




X 


Temp (Matri 


x L f +1) ® 29 






I |2 i 4 


, . X Sj 6. 2, 3,1* 




No records 


§ 29 1 Init. Matrix 






I |2 i 5 


i i i C. Q 1, 2,1 




Matrix Prep c 


tr. J Prep, counter 






1 |2 i 6 


i i i B, Q 0,U,6 







1 






1 |2 i 7 


1 

, , i G| r ■ « > l 




X 




1 J Clear storage area ,,| 






1 |2 i 8 


i i i B | 0, 2 , 7 




Add in matrix 


1 






1 |2 i 9 


, , i a ( q l, i,U 




1 ® 29 








1 I 3 1 o 


i i i Y. Q 0, 2,7 












1 |3 i 1 


i j , E, Q 1, 1,7 




X 


xz 6363 


Mask 


FORM LP-IO 










x 


.0»At N«Btc, J 1 3ft 4 1 X 
CARRIAGE RETURN 
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job no. 


PROGRAM NO. 

F2-129 


PROGRAM PREPARED BY: 

Burggrabe 


PROGRAM CHECKED BY: 

POOL Review 


DATE 

U-18-59 


PROBLEM: 

A Fixed point Matrix and Float 




TRACK 

00 


PROGRAM INPUT CODES 


0. 



h 
l/> 


LOCATION 


INSTRUCTION 


a. 



h 
1/1 


CONTENTS 
OF ADDRESS 


NOTES 


OPERATION 1 ADDRESS 


1 
1 1 1 I 1 1 1 


1 














1 1 1 1 1 1 1 




X 




1 1 1 I 1 1 1 




^J° ,3 ,2 


i , X,S . 6, 2 ,6 p 


/ 


Matrix Lp 


H 




1 1 1 1 1 1 1 






|3 ,3 


i , ,T,0, 0,2 ^ 


/ 








1 1 1 I 1 1 1 






,3 ,4 


i iIiZ , 0,11,0 p 


/ 


BP h - sto 


3 after initialization 


1 ,,..,.. , 






,3 | 5 


, , .R, 0,1,2 U 


' 


X 


•A Matrix P 


rep. 


1 

1 1 1 1 1 1 1 






,3 ,6 


i , ,u , 0,1,* £ 


/ 


Linkage 




1 1 1 1 I 1 1 






,3 ,7 


± l X L Z t A _l ^ 


/ 
1 


— ■> 


BP 8 - sto 


d after matrix prep. 


1 | | | 1 I 1 






I 3 1 8 


8, q 0,T | 0, 2,0 p 


skip matri 


x print out 


1 1 1 | 1 1 1 






,3 ,9 


L , I,B| 6,2,6 1 


' 


X 


(n+1) @ 29 






1 1 1 1 1 1 1 






,4 , 


i i I|C , 6, 2,6 p 


1 


temp 1 - n 


o. of numbers to print 


1 1 1 1 1 1 1 






\* .« 


, | X,B | 6, 2,6 p 


1 




t it ii 


ii ii 


1 1 1 1 1 1 1 






,4 .2 


l l |N| 0,0,5 j2 


1 


1® 25 






1 
1 1 1 | 1 1 1 






,4 ,3 


, il.Aj 6,2,3 £ 


1 


X 


print q 


gives XZ[NNqq] 


1 I 1 I 1 1 1 






|4 | 4 


i i ,1,0,0,^3 


1 




set up Z inst. 


1 1 1 1 1 1 1 






,4 ,5 


, , X,P , 1,6,00 


1 


cr. 






1 1 1 1 1 1 1 






,4 | 6 


1 
, 1 XiZ ■ 0,0,0 p 


1 
1 


delay & ze 


ro (0026 
of no.s 


) 


1 1 1 I 1 1 1 






,4 ,7 


1 
, | X,B | 6, 2,6 ,0 


X 


temp 1 no. 


to print 


,,,;,,, 






|4 | 8 


1 
, , X| S i 6, 2 1 6 1 


1 


(n+1) © 29 












,4 ,9 


1 
i i ,T| 0,0,5 1 


1 


— *■ print tab 






i i i i i i i 






,5 ,0 


. u , 0,0,5 5 


1 




exit tab 1 


oop . 




i i i i i i i 






|5 |1 


, ,I,P, 2,li,0p 


1 


X 


tab 




i i i i i i i 






|5 ,2 


1 
, , X,Z , 0, 0,1 ,6 


1 


delay & 1 


@ 25 I00U2) 


i i i i i i i 






i iS i 3 


, i |A, 0,1,1 k 


1 


1 @ 29 






i i i I i i i 






|5 |4 


1 

i i iU. 0,0 iU 9 


1 
1 

1 




loop 






i i i i i i i 






|5 i 5 


i 
i i i B it i i J 


y 


Lo first n 


o. in row 


i i i i i i i 






|5 ,6 


, , X,R | 0, 6,0 5 










i i i i i i i 






i_ i 5 __lZ 


i xXiUi _Qi_6 ± P£_ 


1 

1 
1 
1 

1 
1 
1 










i i i i i i i 






L_i5 ,8 


r 
_L J_XiZ_|{__L.. J_.j1. 


no & q 


no. of no.s printed 

o. of no. a printed 
Lx prep. 


. i i i i i i i 






l ,5 ,9 


_^^1£+*iU>..&. 


X 


temp 1 


1 h 

i i i i i i i 






1 |6 i 


i i .A| 0,0,5 ..£_._ 


temp 1 r 


i i i i i i i 






1 16 i' 


_L_l.__LY_,._0 1 .0. i 5 ..£.. 


__ i i i i i i i 




_i dL±L 


, , X/B . 6 L 2,6 


i 
i i i i i i i 




i .6 i 3 


i , iUj 0,1,5 ,9 


X 


over matr 


1 1 1 1 _|__l L_ 

FORM LP-12 














X 


CARRIAGE RETl 


JRN 
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JOB NO. 


PROGRAM NO. 

F2-129 


PROGRAM PREPARED BY: 

Burggrabe 


PROGRAM CHECKED BV: 

POOL Review 


DATE \ 

U-18-59 


PROBLEM: 

Given: N records of n variables: prepare A Fixed point matrix & float. 


TRACK 

1 


PROGRAM INPUT CODES 


0. 



t- 


LOCATION 


INSTRUCTION £ 


CONTENTS 
OF ADDRESS 


NOTES 


OPERATION ADDRESS Jj 


1 


/ 
















X 








°l 


1,0 ,0 


, , ,B|[, , ,] ' 


X-i(k+l) @ q 


0^k^n-l 








|0 i 1 


, ■ ,»!,[, , ,] ' 


X.fk+l)^ Q 


0< kin-1 








|0 | 2 


I I 1 U| "li Oi^ ' 


J 










|0 1 3 


1 1 1 1 1 1 1 


:x 




| 








,0 | 4 


1 1 1 1 1 1 1 












,0 | 5 


i . iA|f. i il ' 


k 

7. X.J. @?q 










l° l 6 


, , .Cjh , ,] ' 


Z X.X, @2q 










,0 , 7 


i i i B, Q 1, 0,6 ' 


:x 


^ j 










1° 1 8 


i i i A, Q 1, 2,3 ' 


1 @ 29 










1° 1 9 


i i i Ui Q li 2,6 ' 






,000 1 0,0,1 6 


1 




,1 I 


1,0,0 p , q ( 0,0 i 


1®3 


i| 








,1,1 


i i i | i i i 


■_E 


Data Lo + r 


. + 1 + Data Lj» + 1 








1 1 i 2 


i i i 1 i i i 












1 1 i 3 


I i i 1 i i i 


Data Lo 










il 1 4 


...... ,U 


' 1@ 29 










,1 | 5 


, , , 1 , i ,u 


; 'X 


1® 29 










,1,6 


. i i 1 3 w . w i J 


' XZ 6363 Mask 


ll 








,1,7 


i . i 1 3 w, w i .1 


' 11 n 










,1,8 


i . i . . . . 


' A Matrix Lo 










,1,9 


i i . 1 . i i 


lX 












,2 ,0 


i I i 1 i i i 












l |2 i 1 


i i . 1 i i . 


' -n (Records < 


5 29) 








|2 |2 


1 . f Pi 3» 0, 0,0 


' c.r. (Hex) 










,2 | 3 


...... .u 


•.X 


1 @ 29 & d( 


lay 








,2 | 4 


i i i u i ° °. °i° 


' Exit matrix ] 


irep <r— R here 








,2 i 5 


. 1 1 1 1 1 I 












,2 , 6 


i i T | 9 1 . °|S 












l |2 i 7 


i i i y, q i, o,6 


'X 












|2 , 8 


. i__a B . q 1. 0j0 


' 1 




- 1 








1 |2 i 9 


1 _L_±A(_ 1. Il5 


@ 29 










1 i 3 , 


, , . E | q 1,1,6 


3ww.i 










i |3 i 1 


. , . Yi 1. 0,0 


:X 






FORM LP- 10 












X 


CARRIAGE RETl 


JRN 


DOTAL M««lt. J13&4 1 X 
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JOB NO. 


PROGRAM NO. 

F2-129 


•ROGRAM PREPARED BY: PROGRAM CHECKED BY: 

Burggrabe POOL Review 


DATE 

h-18-09 


PROBLEM: 

A Fixed point matrix and float 


01 


PROGRAM INPUT CODES 


_ 



. 


LOCATION 


INSTRUCTION 


_ 



. 

ID 


CONTENTS 

OF ADDRESS 


NOTES 


OPERATION | ADDRESS 


1 


/ 








1 


data L f +1 


- - 






X 








A 1 .' ,2 


■ , , s |° ililfL 








,3 | 3 


■ , , T , o. 1 . P 


1 
1 


(Back to c 


Dmplete row) 








,3 , 4 


1 1 |B , 0,1,0 I 


>,,.[... 






,3 | 5 


, , |A, 0,1,1 \i 




X 


l @ 29 










|3 , 6 


, , |Y | 0,1,0 I 


1 
1 
1 

1 
1 
1 
1 
1 
1 
1 
1 
1 








1 1 1 1 1 1 > 






,3 , 7 


_L _1_____|. Oj.l.lQQi. . 


_3w_i _. 












,3 , 8 


1 
i 1 iE 1 Oilil 17. 














,3 | 9 


1 
_l J _LS_i_QjJLil-.iL... 


X 


J_ala_L f _±l 
( Begin new 


row' 










I 4 I ° 


1 
. . 1 iT , O1I1O Q.._ 












,4 , t 


1 , iB| 0,1,1 3 


Data Lo 


*\ 










,4 | 2 


x __.jY4_Q_l_O.P_ 






r Initialize 
for new record 


, 1 1 , 1, 1 1 






l 4 1 3 


, , .Yj 0,1,0 JL 


X 












|4 1 4 


j j L_ 4 AlAA_ 

J 1 1 . 1 J L !_._ 














L 4 | 5 




D 








1 4 1 6 


. 1 , B ,to?-i? ,1 


A Matrix L 










,4 1 7 


_____ .1X4 Q_i_>_£ 


X 








1 ::![',! ; 






,4 | 8 


. ■ J| Ol 1 ! , 6 


















,4 | 9 


, , ,B. 0,1,2^ 


1 


-n @ 29 














|5 1 


1 
1 1 1 A | 0, 1,1 $ 


1 
1 


1 @ 29 












,5 | 1 


1 1 1 h | a 1 1 2 ii 


X 














|5 1 2 


, 1 ,T, 0,1,5.4 


1 
1 






J 








1 |5 1 3 


, , 1 U 1 0, 1,22 


matrix pre 


p exit 


-» 








i5 1 4 


j j____Z4______0j0 . 


1 




BP k stop 


afte 


r ei 

«*- 1 


ich record 




■- 




1 1 5 - 1 5 


, , X.R, Q 3,0,8 


1 

1 


X 






J here 






1 |5 1 6 


T 

1 1 X, U | Q 3.0 ,0 














1 |5 |7 


, , l_Ll._O_1i0j0. 


1 
1 
1 

1 
1 
1 

1 














I |5 1 8 


1 

±_ L_ .J j- L._ .1 1... . 














1 |5 1 9 


___. i._... L s,..A.i, 1 ,A 


X 


1 @ 29 


f no .3 to prin X ...... 

out of matrix 
int new line 


'!!',!!!! 






1 16 1 


L i_Xl__i__i 2j6iQ 


jtemp_ _noi c 

1 @ 29 
exit prinl 








I 16 1 1 


1 

_____l_j.__|_Q..__i.1 iL 

8,0,0; f q 2,0,0 

, , ,Uj q o,U,i 


1 1 1 i 1 1 1 






1 16 1 2 








1 16 1 3 


X 


back to p 


f'OHM LP- 12 














X 


CARRIAGE RET 


JRN 
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JOB NO. 


PROGRAM NO. 

F2-129 


PROGRAM PREPARED BY: t 

Burggrabe 


PROGRAM CHECKED BY: 

POOL Review 


DATE 

k-lQ-59 


PROBLEM: 

A Fixed Point Matrix and Float 


TRACK 

02 




PROGRAM INPUT CODES 


0. 




LOCATION 


INSTRUCTION 


a. 

h 
tfl 


CONTENTS 
OF ADDRESS 


NOTES 




OPERATION ADDRESS 




1 


1 


















1 


X 










| 


2 ,o | o 


i , X, Z . 0, U|0 ,0 




BP U stop 












,0 , 1 


i , X , B 6, 2 , 3 ,2 




2q ©29 












|0 ,2 


1 1 iY| Oi 2i0 i9 
















|0 i 3 


i i X|B - 6, 2,3 |6 




X 


A Matrix L 


■ 








1° 1 4 


. , X i A | 6 i 2 i 6 (3 




Tn+l)(n+2) @ 




29 










,0 | 5 


, , l,c ! 6, 2,6 p 




ctr 












,0 ,6 


, , ,B,[ , , , ] 




Z X<X„ @2q 












1° I 7 


i i X,R . 0, 9 ,2 $ 




F 


■*- <j 


I25.0R 










|0 | 8 


, , X,U 0, 9 |0 p 






) 










,0 ,9 


. . .Z,[ , , ,] 




q of no. to fl 


oat. 










,1 I 


, , ,c ,r i i i] 




floated no. s 


bore 












i .i iB 1 Oi 2 i0 6 




X 




w 










i i i A 0, 1,1 k 




1 @ 29 














l i IT | Oi ?i0 6 


















i i |T - 0, 2,1 P 


















i i |E - 0, 1|1 ,7 




X 


3ww,i Mask 


i 


1 








,1,6 


i i X,S | 6, 2,6 p 




Matrix Lp + 1 


1 








,1,7 


. i |T , 0, 2,0 6 
















,1,8 


i , X,Z . 0, 8 ,0 ,0 




BP 8 stop aft 


3r float. 










,1,9 


. . ,u,r . . ■ 1 




X 


Exit sub. 












,2 ,0 


i i i I i i i 




open 












i |2 i 1 






i 












|2 |2 


l i i 1 I i i 
















|2 | 3 


i i i 1 i i i 




X 














,2 | 4 


i i i 1 i i i 










i i 1 4 . 






( 2 i 5 


i i i | i i i 
















,2 | 6 


i i i | i i i 
















|2 i 7 


i i i i i i i 




K_ 






» 


. 






|2 | 8 


i i i 1 i i i 






1 








|2 i 9 


i i i | i i i 


— 














,3 | 


i i i | i i i 














|3 | t 


i i i I i i i 




X 








FORM LP- 10 












X 


.0«al M«tii. J1344I X 
CARRIAGE RETURN 





POOL Program no. - F2-129 
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TITLE: Triangular Matrix Scaling 

AUTHOR ; William F. Burggrabe 

Compumatix, Incorporated 

DATE; September k, 1959 



PURPOSE: Given a matrix in floating point of the form 



N Z 



x i x 9 

—=- Z 2 
C l C 2 


• • • Zi 


*n 


' X l \ 2 , x l 
v C-l ) - C-l 


X ^ 

C 2 


X-j_ X n 
ii • • Zj — — x — — — 

c i c n 



■ (t) 

and scale factors 1, C l3 C 2 ... C n compute the matrix. 



N Z X-j^ Z X 2 EL ... Z X 



n 



X X, Z X-, Xp ... Z X-jX 



2 
Z X„ 



RESTRICTIONS: 



l) The following locations must contain: 

6236 Lo matrix 

6237 Lo floating point 21;. 
62ljl Lo scale factors 

6261 (n+1) @ 29 - matrix order 



Page 1 



POOL Program no. - F2-12? 

2 



2) The scale factors in floating point must be stored 
in the same order as the matrix components. 



3) Normal limitations of 2i|.0, 11.6-12.6. 

U) All scale factors must be stored, (i.e.) where a scale 
factor is unity it must be stored 

X l 
(1 x 1 x N = N, 1 x G. x — ~ = 1 , etc. 

i n-i 1 



c l 



CODING INFORMATION: 



A) Subroutine storage - l|l sectors 

External Storage 6236 Lo A matrix 

6237 Lo floating point 

62J4I Lo scale factors 

6260 Temp storage 

6261 (n+1) @ 29* 

B ) Linkage 

R Lo + 1*0 
U Lo 

C) Input - the following must be stored in memory 

1) A matrix in floating point 

2) Scale factors in floating point 

D) Output scaled A matrix. Stored in the same location as 
unsealed matrix. 

E) Constants 

Lo + 30 XZ6363 mask 
Lo + 31 1 @ 29 

F) Timing (n+l)(n+2) x .65 sec. 



From CI A matrix prep, routines 
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POOL Program no. - F2-129 

2 



Matrix Scaling 



Init. 

'Floating Point 
Lo < Matrix 

Scale Factors 



Enter Floating 
Point 



Compute 



n X,X, 
C xi x C , x Z LJ_ 

X1 X ° 1 C xi C xj 



Inc. Add. 
i'= i + 1 




No 



Yes 



1 



Inc. Add. 
j + 1 = y = i 




No 



Exit 
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JOB NO. 


PROGRAM NO. 

F2-129 


'ROGRAM PREPARED BY: 

Burggrabe 


PROGRAM CHECKED BY: 

POOL Review 


10-U-59 W 


PROBLEM: 

Matrix Scale 


TRACK 

00 


PROGRAM INPUT CODES 


a. 


h 


LOCATION 


INSTRUCTION £ 


CONTENTS 
OF ADDRESS 


NOTES 


OPERATION ADDRESS \j 


1 


1 














1 


X 




1 1 1 1 1 1 1 







°,0 ,0 


i i X,B| 6, 2,3 ,7 ' 


Lo F. P. 




' ' L ' T — ' " 






l0 i t 


i i ,Y 0, 0,1 ,1 ' 












,0 |2 


i i i Y , °i °i 1 i 2 ' 












,0 | 3 


, , X,B, 6, 2,3 |6 ' 


:z 


Lo Matrix 


■ 








,0 | 4 


i ■ |T. 0, 0,1,3 ' 




w 








|0 | 5 


. i iT| Q 0,1 ,7 ' 












,0 ,6 


, , X, B , 6, 2 , k ,1 ' 


Lo Scale Fact 


ors 








,0 ,7 


, i iY, Q 0,1, U ' 


_Z 












,0 ,8 


i i iY, Q 0,1 ,6 ' 












( ,9 


, i XiA| 6 2i6il ' 


(n+l)@ 29 -*• 


Lf +1 (Scale Factors) 








I 1 |° 


, , X,Y, 6, 2,6,0 ' 












,1,1 


. i ,R.[ i , i 1 


~x 


F.P. Linkag 


• 








1 1 i 2 


, i i u | r . i ,i 












1 1 t 3 


i i i pj i i i i 


7. X X 










il i < 


, i iMj[ , , , ] 


C^ J 










i > i 5 


, , x, u, q 0,0,0 


•X 


Ace -^ Mult 










,1,6 


i , ,M,[, , , ] 


, Cx. 


i 








,1,7 


, , ,Hj[ , , , ] 


' Ox -i Qx -i Z X^Ji 


n 








,1,8 


i i X Ei Q 0, 0,0 


' Exit F.P. 










il 1 9 


i . i Bj Q 0,1,3 


^ 


Id of Z X,X 










,2 ,0 


i , , A ( q 0,3,1 


' 1 @ 29 










i |2 i 1 


, , ,T| 0,1,3 












|2 | 2 


i . i Yi Q 0,1,7 












,2 | 3 


, , i Bi 0, lili 


' XI Add. of CV 










,2 , 4 


, , i A| Oi 111 


' 1 (2) 29 










|2 i 5 


i i i Y| Qi l'U 












,2 , 6 


i i , e. q o, 3,0 


' mask 










1 i 2 | 7 


i i X s. 6, 2, 6,0 


lX 


f +1 of C x '. 










1 |2 i 8 


i , i Ti 0. lil 


' loop 


• 








l |2 i 9 


i , i Ui Q 0i 3 i 2 












1 |3 i 


, , X Z, 6. 3, 6,3 


' mask 










i i 3 i 1 


i i X Z. Oi 0,1 


iX 


1 @ 29 




FORM l.P-10 












X 


CARRIAGE RETl 


JRN 


KO'AL H'lll, J 1 H64 1 X 
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PROGRAM NO. 

F2-12? 



PROGRAM PREPARED BY: 

Burggrabe 



PROGRAM CHECKED BY: 

POOL Review 



PROBLEM: 



Matrix Scale 



PAGE 





I 2 



DATE 

10-U-59 



TRACK 

00 



PROGRAM INPUT CODES 



J L 



FORM LP-12 



3 



LOCATION 



0,3 



JJL 



_ll 



Jl 



dl 



Jl 



INSTRUCTION 



OPERATION 



ADDRESS 



B_j_0 



E 



6l 






QJLX 



°Ai_ 



0j3_p_ 



2.6 P 



OjLJL 
kiP_.P_ 



1 



i 1. 



_L 



CONTENTS 
OF ADDRESS 



1 ©29 



loop 



NOTES 



BP U stop a 'te r scali ng 



3 



E 



3 



*K 



/ 



iJ 



CARRIAGE RETURN 
, - ,-.-,Mr..-r.r.MAi ^roB_c^rnL. 
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TITLE: Calculation of Means, Standard Deviations and 
' Correlation Coefficient Triangular Matrix 

AUTHOR : William F. Burggrabe, Jr. 
Compumatix, Incorporated 

DATE: September 23, 1959 



PURPOSE: Given a floating point A matrix as prepared by either the 
~~ fixed or floating point A matrix preparation subroutine, 

compute the means, standard deviation and the all simple 

correlation coefficients. 

n = Number variables 

N = Number records 

The Rij matrix replaces all cross product terms in the 
original matrix. The means and standard deviations (n 
values each) will be stored in the correct order starting 
in the Lo specified. Each value is printed out in normal 
floating point form as computed. (See sample problem) 



RESTRICTIONS : 



1) Normal restrictions of 2J4.O, 11.6-12.6 

2) Required computational storage 
mean Lo thru Lo + (n-l) 

(T x Lo thru Lo + (n-l) 

A Data matrix Lo thru Lo + 



(n+1) (n+2) 
2 



-1 



A Rij matrix replaces all but first row of A data matrix 
3) The following information must be supplied: 

Location 6233 n @ 29 

6236 Lo A data matrix 

6237 Lo floating point 

6238 Lo means 

6239 Lo standard deviations 

6258) 

62^9) Temp, storage area 

6260) 

6261 n+1 @ 29* 

^Supplied by Compumatix, Inc. A Matrix Prep. Progs. 
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k) Output: 2U.0 - 12.6 format 
GENERAL INFORMATION ON METHOD: 





The 


: equati 


ons used are: 






1) 


T i ' 


N 








2) 

i 


Stands 


ird deviation 






J : 


CX 2 
i 


(I) 2 




1 


N 






3) 


Correlation 


coefficient 








2X ± X : 

N 


- -T.T. 

i J 




R U 












CT X . 


^*o 




CODING 


INFORMATION: 









A) Storage - subroutine 2 tracks and 17 sectors 

Calculated storage as described above under 
"restrictions" 2 and 3 

B) No calling sequence required. 

Linkage (R Lo + 216 
(U Lo 

C) Input: None 

D) Output: Format of 21*. - 12.6 

All answers are printed as computed. The correct number 
of tabs preceeds the printing of each diagonal element of 
Rij matrix (See sample problem). The tabs may be eliminated 
by changing t[0209] in 0207 to u[0208] or u0127. 

E) Constants 0021* 1 @ 29 

0201 XZ6363 Mask 

and areas mentioned under Restriction (3). 
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F) Timing: The approximate times including printing are: 
Means: 3 n sec. 

Standard Deviations: I4 n sec. 
Rij matrix: n(n+l) x 2.25 sec. 
Where n is the total number of variables. 

G) Program stops: 

005l Breakpoint 8 T ± calculations complete 
0022 Breakpoint k (T x ± calculation complete 

0215 Breakpoint 8 Rij calculation complete 
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Flow Diagram 
Calculation of Means, Standard Deviations and 
Correlation Coefficient Triangular Matrix 
Floating Point 



0000 



Init. 
Lo F.P. 


1 


' 0007 


Init. 
Lo JMatrix ? 
(& First ZCXj) 


' 


0017 


Init. 
/•Means 
Lo | Standard 
(. Deviations 



* 



Floating Point 
Calculation of 
X, = ZXi & p r int 



Inc. Add. 
i + 1 = i 1 



No 



(bpI 

LaJ 




Yes 



Floating Point 
Calculation 



<xi 



" N 



) ; 



Inc. Add. 
i + 1 = i' 




Yes 



BP 



Floating Point 
Calculation 



ZXjXj 



ij 



X iXj 



xi <J xj 



Inc. Add. 
i + 1 = i 1 





BP 
8 



Yes 



EXIT 



No 



No- 



-*» 



Print 
j' - 1 tabs 
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JOB NO. 


PROGRAM NO. 

F2-129 


PROGRAM PREPARED BY: 

Burggrabe 


PROGRAM CHECKED BY: 

POOL Review 


DATE 1 

U-20-59 


§ 




PROBLEM: - ~~ ~ ~ 

Given: A Matrix; Compute & Print X., , (Tx^ , Rij Matrix (F.P. ) 


TRACK 

00 






PROGRAM INPUT CODES 




h 
1/1 


LOCATION 


INSTRUCTION 


(L 


t- 


CONTENTS 
OF ADDRESS 


NOTES 






OPERATION ADDRESS 






1 




/ 






















/ 


X 














q o l0 ,o 


i , X,B , 6, 2 ,3 i7 




Lo Floating Point 












1 1 t 1 


, , Y 0,0,3 ,3 




mean) , 














1 1 ,2 


i i i Y 0,0,3 k 




mean / 














1 |0 i 3 


i i ,T j 0, o ,5 ^ 




X 0"x, I b 


(1 










1 l" | 4 


. . .Tj 0,0,5 (6 




o-x \ 














1 |0 i s 


i . iY | 0, 1 ,2 ,7 




rij 1 , 














1 i<> |6 


i i ,Y . 0,1,2 J3 




*■ y sub 
rij > 














1 1° 1 1 


i i x i b | 6,2,3 (5 




/\ Lo A Matrix: 












1 1° I 8 


i i iT , 0,0,3 6 




mean ^ 














l 1° 1 9 


i i iY 1 0,0,6 JL 




a 


Lo of N 












1 i 1 | o 


, , , T | o, 1 ^ 3 




ri.i -J 














i l 1 i 1 


i i ,A . 0, ,2 ^ 




X 


1^ 29 


<P 










1 l' |2 


i i iT , 0,0,3 5 




Lo first Z< 


X.) (for mean) 












I 1 1 i 3 


i , XiA | 6,2,3 3 




S ives Lo firsl 


. ^(X. ) 2 (n® 29) 












i il i 4 


i i ,Y 0, ,<S p 




crx sub 


i 












1 |1 i 5 


, i i Y | O, 1 ^ 2 




X 


rij") , 














l i 1 i 6 


i i ,Y 0,1,3 ,7 




c — * SUD 

rij 













1 i 1 i ~> 


i , X,B | 6,2,3 




Lo mean 














1 |< | 8 


i i iY | 0i0i3 ,7 




mean sub 














1 |1 l 9 


, , ,Y , 0. ,5 .7 




X 


CTxilsub 














1 |2 ,0 


, i iY | 0,0,5 e 




crx±) 














l |2 i 1 


i i J 1 0,1,2 p 




rij 








— 1 1 1. 1 III 






1 |2 |2 


, , ,Y | 0,1,3 P 




rij 














1 |2 | 3 


, ,X,P - 1,6,0 D 




X 


cr. 














1 |2 i 4 


, ,XiZ 1 OiOiO I 




delay & 1 <? 


29 












1 |2 | 5 


, iX,B | 6,2 ,3 9 




Lo standar 


i deviations 












1 |2 i 6 


i i iY | 0,1 ,0 




Cxi sub 








— 1 1 1 | 1 1 






l |2 i 7 


, , i Y | o, 1 ^ $ 




X 






9 










1 |2 i 8 


i i ,Y | 0,1,3 6 






( 










1 |2 i 9 


i i X ,B - 6,2 ,3 6 




Lo A Matri 


< 






1 1 1 1 1 1 






1 l 3 l o 


i i X ,A - 6,2 ,6 I 




n + 1@ 29 








1 1 1 1 1 | 






l |3 i 1 


i iX,C , 6,2 ,6 




X 


mean L„ +1 


to test out 




FORM LP-IO 


X 


*ui*i M'aii, J13&4I X 
CARRIAGE RETURN 
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JOB NO. 


PROGRAM NO. PROGRAM PREPARED BY: 

F2-129 Burggrabe 


PROGRAM CHECKED BY: 

POOL Review 


Ji-20-59 




PROBLEM: 

Given: A Matrix- Compute & Print X, . <T^ , Rij Matrix 


00 




PROGRAM INPUT CODES 


a. 



h 

VI 


LOCATION 


INSTRUCTION q 


CONTENTS 
OF ADDRESS 


NOTES 




OPERATION | ADDRESS £ 




1 




/ 




















X 












OJ 


3 ,3 ,2 


, ,X,P - 1,6,0 p ' 


cr. 














,3 ,3 


, , .R ,[ i I I 1 ' 


| F. Point 














,3 ,4 


, , iu ,E i i i 1 ' 


J 


__ 


• 










,3 | 5 


, , >B j[ i i i ] ' 


X_ 


2 X, 












,3 ,6 


i , ,D ,[ , , i ] ' 


No. of recDrds 












,3 ,7 


i , ,H ,[,,,] ' 


- Si 






1 1 1 1 1 1 








,3 ,8 


, ,XiP | OiOiO o ' 


Print X., 














|3 |9 


1 1X1E | 0,0,0 p ' 


_K 


Exit floa :ine point ... .. 












|4 | 


i i iBC 0. i3 7 ' 


Add I ± 














\* l> 


, , ,A | 0,0,2 U ' 


1 @ 29 














,4 ,2 


, , ,Y , 0,0,3 7 ' 






m 










,4 ,3 


, , ,BJ 0,0,3? ' 


_x 


Add ZX t 




w 










|4 l 4 


, i ,A, 0,0,2 k 


1 @ 29 














|4 | 5 


, , J, 0,0,3 $ ' 
















,4 | 6 


i i iE , 0,2,0 I 


3ww.i 


mask 












|4 ,7 


, , XiS | 6i 2i6 


' 7\ 






• 


,111. 








|4 | 8 


. ■ iT, 0,0,3 3 


loop back 


i times 


1 1 1 1 J u 








,4 ,9 


« , X,P | 1, 6,0 ,0 


' cr. 














i ,5,0 


i i X,Z | 0, 0,0 ,0 


' delay 














i i5 , 1 


, , X,Z 0, 8,0 ,0 


IX 


BP 8 stop 


after means 












i 1 5 ,2 


, , X,B . 6, 2,6 p. 


' (n+l)@ 29 














i |5 i 3 


, ,1,0, 6,2,6,0 


i Temp 1 (Ct 


p) 












1 |5 |4 


, , X,P I 1, 6,0 ,0 


' A Matrix I 


r + l(Set for <Tx t ) (Temp : 












1 |5 | 5 


i i i R |L i i i J 


lX 


F.P. 














i )5 i 6 


, , ,u,[ , , i ] 


/ 














l |5 i 7 


, , iPjT i i i 1 


f i 














i 15,8 


. i iM |T i i i 1 


' ^ _ 














l |5 i 9 


i , X, H . 6, 2,5 ,9 


lZ 


Temp 2JX, 


I 2 __. - . 


• 










1 |6 i 


1 1 1 B 1 [ 1 1 1 J 


i 


n ~2 

z x i 




1 1 1 1 1 1 








A^^ 


, , ,D|[ i . , 3 


No of rec< 


Drds 












j |.6 1.2 


, ,1,8,0,0 9 q 


' M ._ 




1 L_.J_ f- 1 L 








i |6 i 3 


■> ; 


V~~ 






1_L_J L_4_J L 

FORM LP-12 












A 


•> 

CARRIAGE RET 


UHN 
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JOB NO. 


PROGRAM NO. 

F2-129 


PROGRAM PREPARED BY: IPROGRAM CHECKED BY: 

Burggrabe POOL Review 


DATE L 

U-20-59 '1 


PROBLEM: 

Given: A Matrix; Compute & Print 1., <r x± , Rij matrix (F. P. ) 


TRACK 

01 


PROGRAM INPUT CODES 


tt. 



h 
11) 


LOCATION 


INSTRUCTION 


a. 



h 


CONTENTS 
OF ADDRESS 


NOTES 




OPERATION ADDRESS 


1 


1 


















X 








J1,0 ,0 


■ •■ ,h ,r , , , i 




Standard Deviation 








,0 , t 


i ,X,P 0, ,0 p 




Print Stan. 


Dev. 










,0 | 2 


i 1X1E | OiOiO 














,0 | 3 


i i iB | 0,0,5 ,7 




X 


Add. of X 4 




-II 








,0 | 4 


i i i A °i°i 2 h 




1 @ 29 










|0 | S 


i i J °i °i^ i 7 














,0 ,6 


i i ,T . 0. 0,5 J3 














|0 | 7 


i i |B | 0,1,0 p 




X 


Add. ofCfXi 




1 1 1 I 1 1 1 






,0 ,8 


i i i A ! ! i 2 h 




1 § 29 










,0 | 9 


, , ,Y | 0,1,0 p 














,1 ,0 


, , X,B - 6,2,6 P 




Temp 1 (n+1 


,n,n-l,etc. ) 








,1,1 


i i i s 1 0, ,2 fr 




X 


1 O 29 


1 








|1 1 2 


, , X,H 6, 2,6 ,0 




Temp 1 










1 1 i 3 


i i i A . 0, , 6 ,0 




n 
Add. of ZX, 


<L 








il l 4 


i i |T | 0, 0,6 ,0 














ll i 5 


i i i B | 0, 1 2 ,li 




X 


1 @ 29 










, 1 ,6 


i i XiS | 6i 2,6 ,0 




Temp 1 (Tes 


t out when = l) 1 








|1 |7 


i i i T 1 °i °i£ 5 




loop back ( 


a times) 








,1 i 8 


, , X,P, 1, 6,0 ,0 




cr. 










,1 , 9 


, | X,Z 0, 0,0 ,0 




XL 












,2 ,0 


, , X,P 1, 6,0 ,0 




cr. 










i |2 i 1 


, | X,Z . 0, 0,0 ,0 














|2 | 2 


i i XiZ | Q liiO iO 




BP k stop a 


fterCx, scale 








i2 |3 


, , X,B. 6, 2,3 ,8 




X 


Lo mean 










|2 i 4 


i i X.Aj 6, 2,3 .3 




n @ 29 










|2 | 5 


, , X, C . 6, 2,6 ,0 




mean L„ + 1 


(Temp 1) 










|2 | 6 


i i I|C| 6, 2,5 ,9 




Tab counter 


(Temp 2) 








|2 i 7 


, , ,R|[ , , J 




X 


F.P. 










l2 i 8 


. , .Ujt , , ,] 






< 








1 1 2 ,9 


1 , .Pit 1 1 ll 




X. 










,3,0 


, , ,M,r , , ,i 




1 
X. 










l |3 i 1 


i i X, H - 6 2,5 ,8 




X 


Temp 3 


(xX) 


FORM LP-tO 














x 


CARRIAGE RETL 


x J 

IRN 




JI3&4IX 
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JOB NO. 


PROGRAM NO. 

F2-129 


'ROGRAM PREPARED BV: PROGRAM CHECKED BY: 

Burggrabe POOL Review 


lt-20-59 


PROBLEM: 

Given: A Matrix ; Compute & 


Print X i3 <Tx± , Rij Matrix (F.P.) 


01 


PROGRAM INPUT CODES 


a. 


t- 


LOCATION 


INSTRUCTION 


0. 




CONTENTS 
OF ADDRESS 




NOTES 


OPERATION | ADDRESS 


1" 


/ 








/ 














X 













A 1 ^ .2 


,', ,B,[, , , 1 


Z X-jXj 








,3 | 3 


■ , ,v. ■ - ] 


/ 
1 





No of reccrds 











,3 ,4 


i i X,S 6, 2,5,8 


Temp 3 




| 






,3 | 5 


i , ,D [, , , ] 




X 


<Tx . 




' 






,3 | 6 


, , ,y, , , 1 


/ 
/ 




5"x . 

_ i 










,3 | 7 


i . . H, [ , , , 1 


rij 












,3 ,8 


i i X P| OiOiO 


Print ri; 










,3 | 9 


i i Xi E. Q 0, 0,0 


Xl Exit F.P. 












,4 | 


i i , B, q 1,3,2 


/ 


add Z X.,X. 












,4 | 1 


i i iA| q 0,2 ,1| 


/ 
/ 
/ 
/ 
1 




l@ 29 










I 4 1 2 


, , , T| q 1,3.2 






►,,,,,, 






,4 , 3 


, , iTj q 1,3,7 


X 














|4 | 4 


i i i B| Q Ii3.6 




add tTx-i 










(4 | 5 


I i I Ai Oi ? i)| 


-JL®. 29. 










,4 i 6 


i i i Yj q 1,3 ,6 


1 
1 
1 














,4 | 7 


, , ,B| Q 1,3,0 


X 


add JL 




i : : : ; : : : 






,4 | 8 


i i i A | 0. 0.2 i)| 




1 @ ;>q 










|4 i 9 


, i iT| 1, 3 iO 


1 












,5 L 


i , i F,| 2i Oil 


1 
1 
1 


3wwj mast 












,5 | 1 


, i X S| 6 2, 6,0 


X 


mean Lp +] 












,5 i 2 


, i , T| q 1,2,7 












i |5 i 3 


, , ,'B| q 1,3,5 


1 


add <^'x i 










,5 | 4 


i i i A | Oi 2 ill 


1 
1 




1 @ 29 




- -. - 








|5 | 5 


■ . i T| q 1,3,5 


X 












i |5 i 6 


i i i Y| n ii 3 ifi 


1 












|5 i 7 


, , ,b, q 1,2,9 


1 
1 
1 

1 
1 
1 

1 


add X. 










| |5 i 8 


i i i a, q o, 2 ,U 


~i 
1 @29 










1 |5 i 9 


. i i Y| 1.2 ,9 


X 








i> , , , 1 , , , 






1 |6 i 


, , >Y| q 1.3..1P. 


3wwj mas' 
mean L*. 


+1 


■ - 


i i i i i i i 






i |6 i 1 


. . i E, q 2,0,1 








l |6 i 2 


i , Xi S| 6 2i 6i0 








i |6 i 3 


, i ■ Tj q 2,0,2 


\ / 






FORM LP- 12 














X 


CARRIAGE RET 


JRN 


«, ■■i m i,i . ji :,s4 1 x 
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F2-129 


PROGRAM PREPARED BY: 

Burggrabe 


PROGRAM CHECKED BY: 
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DATE 

U-20-59 


PROBLEM: 

Given: A Matrix; Compute & Print X.. , (Tx± , Rij Matrix (F. P.) 


TRACK 

02 


PROGRAM INPUT CODES 


0. 


LOCATION 


INSTRUCTION 


Q. 


h 


CONTENTS 
OF ADDRESS 


NOTES 


OPERATION ADDRESS 


1 


/ 
















/ 


X 








q 2 l0 10 


1 1 iU| Q 2,1,3 




Exit Path 


-*- 






1 1 1 1 


1 , X,Z 6, 3|6,3 




Mask for ey, 


tract 's 






I |0 |2 


, 1 X,P L, 6,0,0 




cr. 








1 |0 1 3 


1 1 x i B , § 2 \$ |9 




X 


tab ctr. 


1 ! 






1 1° |4 


1 1 1 A - Q 0,2 ,ii 




1 @ 29 


(increase tab ctr 






I |0 |5 


t 1 X |C, 6, 2,5 ,9 












1 |0 ,6 


. 1 XiS| 6 2i* .9 




make ace. n 


egative 






1 |0 ,7 


. 1 1 T , Q 2 , ,9 




X 


Print tab 








1 1° 1 8 


1 1 , U . 0, X , 2 ,7 




Exit tab lc 


op 






1 1 t 9 


1 , X,Z | 0, 0,0 ,0 




delay 








1 I 1 |° 


, , XiP, 2, lnO ,0 




tab 








1 1 1 t ' 


1 .1 1 A | 0, , 2 h 




X 


1 @29 








I 1 1 1 2 


1 1 , U | 0, 2 , ,7 




loop 








1 1 1 1 3 


1 , XiP. 1, 6,0 ,0 




cr. 








1 i1 1 4 


1 , X|Z, 0, 0,0 ,0 












I |1 1 5 


, , X,Z | 0, 8,0 ,0 




X 




i 






1 (1 1 6 


1 , iU|f ,,,1 






\ 






1 i< 1 7 














1 ,1 1 8 














I (1 1 9 






X 










1 |2 , 














1 (2,1 














1 |2 |2 














I |2 1 3 






X 










1 |2 1 4 










,111,1! 




1 (2 1 5 














1 |2 , 6 














1 |2 1 7 






X 










1 |2 1 8 














1 (2 1 9 














1 |3 1 










1 1 1 | 1 1 1 




1 |3 1 1 






X 




















V 
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TITLE: Convert Triangular Rij matrix to a square Rij matrix - modified. 



AUTHOR : William F. Burggrabe, Jr. 
Compumatix, Incorporated 



DATE: 



October 13, 1959 



PURPOSE : Given a triangular matrix as prepared by the subroutine 
'"Calculation of Means, Standard Deviations and A Correla- 
tion Coefficient Matrix", convert the A correlation 
coefficient portion to a square Rij matrix with the 
initial location of A Rij equal to the initial location 
of the D Rij. This routine also sets automatically 
location 621+0 to the correct address. 



RESTRICTIONS : 

A) The Lo A Rij matrix = Lo D Ri j matrix 

B) The entire subroutine is in fixed point and will 
operate on both fixed and floating point format 
matricies, but will not operate on matricies in 
an extended range format. 

C) External storage as on attached sheet. 



METHOD: The subroutine takes advantage of the symmetry of the 
matrix and sets cell ij = cell ji. A looping process 
builds both the row and column cells of the square matrix 
while stepping through the matrix from the bottom. 



CODING INFORMATION : 

A) Storage: 56 sectors 

B) Constants: Lo + J?l 1 @2 

Lo + 52 1@29 
Lo + 53 Mask 

C) Linkage: R [Lo + ^] 

U [Lo ] 

D) Stops: Break point k at end of conversion 

E) Timing: Approximately n(n+l) seconds 

6 
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PROGRAM NO. 
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PROGRAM PREPARED BY: (PROGRAM CHECKED BY: 
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DATE 1 

10-10-59 


PROBLEM: 

External Storage A — *• n Rij Modified 


TRACK 

62 


PROGRAM INPUT COOES 




i- 


LOCATION 


INSTRUCTION 


0. 



h 


CONTENTS 
OF ADDRESS 


NOTES 


OPERATION | ADDRESS 


1 




/ 


















/ 


X 










6 , 2 ,3 ,2 
















1 |3 |3 


, ,x,z ,r .ill 




n © 29 


No variables -a- 








, |3 ,4 


















1 |3 |5 






X~ 


II 








1 ,3 ,6 


, ,I|Z,[ , , , ] 




Lo A data Matrix * 








1 ,3 ,7 
















1 |3 i 8 














1 I 3 I 9 






X 










1 1 4 | 


i |X,Z,[ , , ,1 




Lo □ rij ma 


brlx 









1 l< l< 
















1 l< |2 
















i ,4 ,3 






X 




1 








1 M | 4 






Note: -tf Thes 


3 locations 








l |4 ,5 






must be fille 


i prior to 








i |4 i 6 






entry into th 


3 subroutine 








1 |4 ,7 


1 1 1 1 1 1 1 




x 












1 |4 | 8 


1 






II 








l ,4 ,9 






Note : / The 


se locations 








, ,5 ,0 






normally fill 


sd by A Matrix 








1.5,1 






X 


Prep . sub . 










, ,5 ,2 
















l |5 i 3 
















1 i5 |4 
















1 iS IS 






X 












l iS ,6 
















1 iS |7 
















1 i5 18 
















i iS .9 




X 


Temp, s 


iorage > 








1 |6 i 








( 








1 16 ,1 






n + 1@ 29 


.* y 








1 l« |2 
















, ,6 ,3 




X 


(.n+lKn+^J 


® 29 * J 


? 


FORM LP- 12 










1 


X 


M:...l M Imi Jl 144 1 X 

CARRIAGE RETURN 
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F2-12° Burggrabe 
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PROBLEM: 

^ ri i „^ D v i^ Matrix Conversion 


00 


PROGRAM INPUT CODES 



h 
VI 


LOCATION 


INSTRUCTION q 


CONTENTS 
OF ADDRESS 


NOTES 


OPERATION ADDRESS Jj 


j, 0,0 ,0 ' Lo , , 


/ 












/i °i ° i° i L ? i i 




X 








q o,o ,o 


i , X,B, 6, 2,3,6 ' 


Lo A data natrix 


— i — i — i — | — 1 — I — I — 




1 i0 , 1 


, XA 6, 2,6,3 / 


(n+l)(n+2j,, 0Q 

<-> — 








,0 ,2 


, , ,S' 0,0,52 , 


1 @ 29 




.,.[... 






,0 ,3 


j\ 0,1,7 ' 


X 


B[ ] 


Lf Adata matrix 








|0 i 4 


IB 1 6236 , 
i i i i i i i 


Lo A data natrix 








,0 | 5 


X a' 6 2 6 1 / 
i i i i i i i 


n + 1 @ 29 










,0 ,6 


, , X ,C 6, 2,U,0 ' 


Lo A-H i matri: 


& Lo O ri,1 matrix 








,0 ,7 


IC 62 60 / 
i i i i i i i 


1 


Temp 1 -*0 


*v 








|0 I 8 


. ■ X , C | W? ' 


Temp 2^0 


Init. ctr's 








,0 ,9 


i , x , c! 6, 2,5,8 ' 


Temp 3 -* 


i 








,1 ,0 


, I X , B I 6 , 2 , 3 , 3 ' 


_. „ n ® 29 


; 


! i i i | i i i 






,1 ,1 


I I X I N I 6 . 2 . 3 I 3 J 


_E 


n@ 29 


I 








|1 |2 


, i , M '| q o,^, 1 




1 @ 2 — ^ 


n @ 29 ! 








|1 |3 


i , x i A , 3 2 , ^i° 


Lo D rij mat: 


•ix 








,1 |4 


, , .Sj q 0,5,2 


' 1 © 29 -» 


Lf., n rij matrix 








,1 ,5 


, i i T. q 0,1,8 


X 


H [ ] 




i :::;:: : 






,1 ,6 


Y 19 

i i i I i i i 


c [ ] 










J .7 


B [ ] 
i i i | u i i i 


' ri.i A matr 


.X 








1 .1_jJ 


h'[ ] 

i i i ii i i 


' rij n (row 


) 








l il |9 


, , ,cj[, , , ]■ 


lX 


r.ii n (col 


) 








l |2 i 


, , ,B. Q 0,1,7 


B [ 1 










i |2 il 


, . , Sj 0,5,2 


' 1 @ 29 










1 |2 |2 


, , il: 0,1,7 


B [ ] 










l .2 |3 


B 18 
i i i 1 i i i 


lX 


H [ ] 










l |2 i 4 


,,, s ; qo,*, 2 


' 1 @ 29 










i |2 |5 


,t| q 0,1,8 


H [ ] 










i |2 , 6 


, , XT, 6, 2,5,8 


' Temp 3 










i iZ .7 


, , , B, Q 0, 1,9 


lX 


c [ .] 




>!!!'!!! 






1 |2 i 8 


, , x, s ( 6, 2,3,3 


i n® 29 










i |2 ,9 


, , ,T. Q 0,1,9 


c [ ] 










l .3 ,0 


, , , E | 9°,5,3 


' 3wwj 


mask 








i i3 ,1 


, , JjSJ 6 2,5,8 


^ 


« Temp 3 




FORM LP-IO 












X 


CARRIAGE RET 


URN 
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PROBLEM: 

A rij —> D rij Matrix Conversion 


TRACK 

00 


PROGRAM INPUT CODES 


0. 



h 
VI 


LOCATION 


INSTRUCTION 


0. 



h 
V) 


CONTENTS 
OF ADDRESS 


NOTES 


OPERATION 1 ADDRESS 


1 


1 
















1 


X 












0, 0,3 ,2 


i i , T | 0, 0,3 k 




test out 










1 |3 i 3 


i , r U 0,0,1 ,7 




loop 








l |3 | 4 


! 1 | B i°| 1 ? 




c[ ] 








l |3 i 5 


i ,X,A 6,2,5 ? 




X 


temp 2 


0,l,2,...n-l J 






1 i 3 , e 


i i J °i Oi 1 P 




H[ ] 


1 






1 I 3 i 7 


i iX,B | 6, 2 6 9 




Temp 2 


r i 






1 l 3 i 8 


i i iA | 0,0,5 ,2 




1@ 29 


/inc. ctr 2 






1 1 3 | 9 


i i ZiC | 6i 2,5 p 




X 


temp 2 


. 






I |4 , 


i i XiB ! 6i2i5 6 




= old HT 1 








1 I 4 1 ' 


i i X,A | 6, 2,6 P 




temp 1 


0,n,2n,...(n-l)n 






j ,4 j 2 


i i iY | OiOil 9 




or i 








i l 4 i 3 


i . X.B - 6, 2,6 ,0 




X 


temp 1 


1 






1 |4 i 4 


1 
i i X,A, 6,2,3 3 




n @ 29 


}inc. ctr 1 






1 |4 i 5 


"I 

i i XiC | 6i 2,6 




temp 1 


f - 






i |4 i 6 


i i iB | 0, 0,1 ,8 




H[ 1 








1 |4 i 7 


i i |E | 0, 0,5 3 




X 


3wwj 


mask 






l |4 i 8 


, , X,S 6, 2,i| ,0 




Lo n rij m 


atrix 1 






l |4 i 9 


i i iT j 0, OS h 




test out 








1 |5 i 


i i iU , 0,0,1 ,7 




loop — > 




,000 ( § ^ 


1 


1 |5 i 1 


2 0, 0,0 , 0, 0,0 p 




X 


1 © 2 








l |5 i 2 


, , , ! , , h 




1 & 29 








i |5 i 3 


i i i 1 3i w,w i.i 




mask 








1 [5 | 4 


i i ?z j o, U,o p 




BP k sto ( > 


i!'t,(>r conversion 






1 |5 | 5 


i , ,u,[ , , , 1 




X 


exit 




. ' 


1 


l |5 i 6 














1 |5 i 7 














l |5 i 8 














1 |5 i 9 






X 










1 |6 i 


i i i 1 i i i 






( 






l |6 i 1 














. .6,2 














_L |6 i 3 






X 






FORM LP-12 














X 


CARRIAGE RETL 


RN 


NCfAt H-WI I JIJ14 1 X 
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TITLE: Compute and Print Matrix Inverse 

AUTHOR : William F. Burggrabe, Jr. 
Compumatix, Incorporated 

DATE : April 21; , 1959 

PURPOSE : A) Initialize and transfer control to the matrix 
_____ inversion routine (29.0). 

B) Printout the inverse if it is desired (Transfer 
control up) 

RESTRICTIONS : 

1) Normal restrictions of 21;. 0, 12.6 and 29.0 

2) See attached sheet for external information that 
must be supplied. 

3) Output - 12.6 format 
CODING INFORMATION : 

A) Storage: 55 sectors and under two above. 

B) Linkage: R (Lo + 5U) 

U (Lo ) 

C) Input: None 

D) Output: 12.6 format with a carriage return after 
each row. 

NOTE: Transfer control down skips printing to Lo + 53. 

E) Constants: 



Lo + 28 
Lo + 29 
Lo + 30 
Lo + 31 



1 @29 
lli @ 29 
XZ01U9 

l @ lU 



F) Timing: 

Inversion: Approximately 1.08 ir seconds 
Printout: Approximately 2 n seconds 

G) Program Stops: 

Lo + 15* B.P. 8 inversion complete 
Lo + 53 B.P. I4 printing complete 

^Depression of transfer control before continuing causes a transfer to 
Lo + S3>, thus eliminating the printing phase of the program. 
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PROGRAM CHECKEO BY: 

POOL Review 


DATE ^^H 
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External Storage 


TRACK 

62 






PROGRAM INPUT CODES 


Q. 

1- 
1/1 


LOCATION 


INSTRUCTION 


B. 



h 
in 


CONTENTS 
OF ADDRESS 


NOTES 






OPERATION ADDRESS 






1 




/ 




















/ 


X 














6,2,3 ,2 





















1 |3 ,3 


, .x.-z^ , , , ] 




n@ 29 = Matrix Order # 










1 |3 ,4 




















. ,3 ,5 






x - " 


1 


■ 










1 |3 |6 




















. .3 ,7 


1 1 x,z,r , , , 1 




Lo floatin 


g Point -if 













I i 3 1 8 


















1 ,3 ,9 


1 1 1 1 1 1 1 




xr ' 














I I 4 1 ° 


, , x.z'c , , , ] 




Lo D rij 


matrix % 












1 |4 ,1 


1 1 1 1 1 1 1 


















1 I 4 .2 


, , X.ZJL , , , ] 




Lo inversi 


on sub. 29.0 * 












,,4,3 






X 


— - 


~ ip 










1 |4 1 4 




















1 |4 ,5 






* Note: These locations 












1 ,4 | 6 






must be filled prior 












1 ,4 ,7 






X 


to entrv into subroutine 












1 ,4 , 8 








. ( 


I 










1 ,4 ,9 




















1 ,5 ,0 




















1 ,5 ,1 






X 
















1 ,5 ,2 




















1 |5 | 3 




















1 ,5 ,4 


1 1 1 1 1 


















1 ,5,5 






x 
















1 ,5 ,6 




















1 ,5 ,7 




















1 (5 , 8 




















1 |5 ,9 


r. ...... 1 




X 


Temp & ctr 




• 










1 |6 | 


[.,,,,,] 




Temp & ctr 






1 1 1 1 1 1 1 






1 ,6 ,1 




















1 ,6 ,2 














— 1 1 1 1 1 1 1 






1 ,6 ,3 






>^ 








FOftM cs* :: 


X 


CARRIAGE RETURN 
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DATE 

U-2U-59 


PROBLEM: 

Call-in, Compute & Print - Inverse (Rice F. P.M.I.) 


00 


PROGRAM INPUT CODES 


0. 



t- 


LOCATION 


INSTRUCTION 


a. 


h 


CONTENTS 
OF ADDRESS 


NOTES 


OPERATION ADDRESS 


; ' Lo 


/ 














/ Lo 




X 








0_j 


,0 ,0 


, , X i B | 6 i 2 ,3 r 3 




n@ 29 










|0 i 1 


. i .D 0,0,3 JL 




1 @ iU 










,0 ,2 


i , X|A . 6| 2 iU P 




Lo n ri.i mati 


ix 








i0 | 3 


i i iC | 0,0,1 k 




X 












,0 | 4 


, ,X,B , 6,2 ,U ,2 




Lo matrix Inversion sub. (Rice M.I.) 








,0 ,5 


i . iY 1 OiOil 3 














,0 ,6 


i i |A | 0,0,2 9 




XZ 1 k 










,0 | 7 


, i ,Y | 0,0,1 2 




X 












|0 ,8 


i i iA | 0i0i3 




XZ 1 k 9 










,0 ,9 


i i iY 1 0,0,1 1 














,1 I 


, ,X,B | 6,2,3 i7 




Lo Floating 


Point 








,1,1 


i.i iY|l" i . .1 




X 


(Matrix Inv« 


rsion Lo + 163) 








1 1 i 2 


, , ,R [ , , ,1 




V Matrix li 


version Linkage 








1 1 i 3 


, , ,ur, , ,] 




J 










il l 4 


[,,,;,,,] 




Code word 










,1 ,5 


i i X, Z, Q 8,0 ,0 




X 


BP 8 stop < 


.fter inversion 








,1 |6 


1 1 X I C , 6 I 2 I 6 I° 




Clear ace. 










|1 i 7 


T q o,t o, 0,5 ,3 




Skip inverse 


print -out j 


" — *- ' — - 1 — 1 — ' — ' — ' — 






,1,8 


• i I|C, 6, 2,6,0 




ctr 1 -*- 


| 








,1,9 


, , X, C ' ^ 2,5 ,9 




XI ctr 2 -> 


i 

i 








,2 ,0 


i i X, B 6, 2,3,7 




Lo Floating ] 


'oint J 








i |2 i 1 


, i iTj Q 0,3.2 






V Init. Print 








|2 |2 


i i . Y | 9 0, 3,3 






( 








,2 | 3 


, , X,B, 6, 2, U ,0 




X. 


Lo Rij M 


itrix Inverse 








l |2 i 4 


, , t y q o, 3 ,U 






i 








|2 | 5 


, , X| Pi 1 6,0,0 




C.R. 


j 








1 |2 i 6 


, , i, z. q 0,0,0 




delay 










l |2 ,7 


, i , tj. q 0,3,2 




X 


transfer o 


rer constants 








1 |2 i 8 


i i x, z | q 0,0,1 




1 @ 29 










l |2 |9 


, i Xi z, q o,i iU 




Ilia 29 










1 |3 i 


. . x, Z| q i,U,9 




Used to set u 


(Lo + lo3 in 
5 ]H1cp F.P.M.T. 








i |3 i 1 


, , x, y. q 0,0,0 




X 


1 @ 111 


1 


FORM UP-IO 














X 


CARRIAGE RET( 


JRN 


huial M«tie, JI3641 X 
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PROBLEM: 

Call-in, Compute & Print - Inverse (Rice F.P.M.I.) 


TRACK 

00 




PROGRAM INPUT CODES 


0. 


t- 
(/> 


LOCATION 


INSTRUCTION 


n 


h 


CONTENTS 
OF ADDRESS 


NOTES 




OPERATION j ADDRESS 




1 


1 












ge 






1 


X 












|3 i 2 


, i ,R|[ ,,,] 




F.P. Linka 










|3 | 3 


, , ,u,[ , , , ] 










1 1 1 1 1 1 1 






|3 | 4 


, , ,B,[ , , , ] 




No. to pri 


nt 










,3 | 5 


, | X, P 0, , ,0 




X 


Print No. 


j 


| 








|3 | 6 


i , X, E | 0, 1 ,0 




Exit F.P. 












,3 | 7 


i i |B, 0, 0,3 h 




add. of nc 


. printed 










,3 , 8 


i i i A 1 Q 1 2 18 


1 @ 29 












,3 | 9 


i i ,Y| Q 0,3 ili 




3 














,4 | 


, , X, B 6, 2 , 6 ,0 




ctr 1 












,4 ! 1 


i i , A ■ 0, , 2 ,8 




1@ 29 












,4 , 2 


, , X,H. 6, 2,6 ,0 


— - 


ctr 1 












,4 , 3 


i , l| S - 6, 2 , 3 ,3 


k 


n@29 


i 


1 








|4 | 4 


i i i T , Q 0,3 ,2 


loop n tim 


es 










|4 | 5 


I , X,P, 1, 6,0 ,0 




cr. 












,4 | 6 


i i x i c | 3 2, 6 ,0 




ctr l-*0 












,4 i 7 


i i X,B| 6 2|5,9 




X 


ctr 2 












,4 | 8 


1 A ! | 1 2 I 8 




1@ 29 


1 


1 








,4 | 9 


1 1 X, Hi 6 2,5 d 




ctr 2 












,5 , 


, , X, Z. Q 0,0,0 




delay 












|5 i 1 


, i I.S, 6 2,3,3 




X 


n@ 29 










1 


,5 , 2 


i i i T , °i °i- i 2 




loop n tiir 


es 










15 i 3 


i i X, Z I Q Ui ,0 




BP U stop 


after printing 










|5 | 4 


i , ,u |L , , , ] 




exit 












|5 I 5 






X 














|5 i 6 


















|5 | 7 


















|5 ,8 


















|5 i 9 






X 














I 6 i 








4 


1 








16 | 1 


















,6 | 2 


















,6 , 3 






X 








FORM LP-12 














X 
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TITLE ; Calculation of Beta Weights, Regression Coefficients ; Partial 
Correlation Coefficient and Standard Error of the Independent 
Variables; the Constant Term (b ) and its Standard Error;the 
Sample Multiple Correlation Coefficient and Standard Error of 
Estimate; the Universe Multiple Correlation Coefficient and 
Standard Error of Estimate. 



AUTHOR : William F. Burggrabe, Jr. 
Compumatix, Incorporated 

DATE : May 26, I960. 



PURPOSE : Given the inverse of the correlation coefficient matrix, the 
means and standard deviations of the variables; compute and 
printout the above. 



RESTRICTIONS : 

1) The last row and column of the matrix contain elements 
corresponding to the dependent variable. (Explanation 
under "Method" ) . 

2) Normal restrictions of 2J4.O and 12.6. 

3) See attached sheet for external storage of information 
that must be supplied. 

h) 25. OR in 0900. 



METHOD: 



A) Given a correlation coefficient matrix of the form: 





X l 


x 2 


x 3 


• • • 


Y 


\ 


R ll 


R 12 


R 13 


• • • 


R iy 


h 


R 21 


R 22 


R 23 


• • • - 


R 2y 


h 


R 
31 


R 
32 


R 
33 


• • • 


R 

3y 


• 
• 


• 
• 
• 


• 
• 


• 
• 




• 
• 
• 


Y 


Ryl 


R o 






V 
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An inverse of the (Rij) matrix is computed yielding: 



V l»23 ... y 

Pl2.3lt ... 
V l,23 ... y 



P21 t 3lj ... 
v 2,13 ... y 



V 2,13 ... y 



. Ply, 23 ... P2y,13 ... 



Pyl.23 .. 
v y,l23 ... 

Py2,l3 .. 



v 



y,i23 ... 



^1,23 ... y V 2,13 ... y V y,123 ... 

B) Where: 

1. The p's, Beta weights, are equal to the coefficients of the equation: 

IL± = H %*2l + P 2 i™h + ... + Pn-1 (*-*)n -l 

^xn-l 



*~xl 



0~x2 



Y' = Predicted Value 

P2 =_Py2,13 ... 
xi, y = means 

^xi'^v = s " tandard deviations 



The p's are obtained by dividing each column element by the 
principle diagonal element in that column. 

NOTE: All p's for all multiple regression equations 

(i.e., 

are available in the inverse. 
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2. The regression coefficients are obtained from the Beta weights 
as follows: 

b = Y - X x £L B x - L $J- P 2 - *3 ^L ^3 ••• 



"2 — p 2 " A 3 
<T X 1 <T x2 J 0"x3 



b l = Pi 



0% 



^xl 



b 2 = B 2 y > etc - 
<5~x2 

3. The partial correlation coefficients are obtained by: 

/ Pab,cde ... I ^ / Pba,cde 
\ "a,bcde ... / \ b,acde 

ab , cde . . . '. * J — 



r 



V. 



V> 



'a,bcde ... / \ b,acde 
However, because of the symmetry of the inverse, 

"ab,cde ... _ "ba,cde ... 
v a,bcde ... ^b,acde ... 



There, each partial correlation coefficient may be obtained by 
dividing the corresponding off principal diagonal element by 
the square root of the product of the row, column elements of 
the principal diagonal. 



I;. Standard error of regression coefficient — 



(l-rjy,123 ...) 
S bi = S • 'V N<r ^ ^- R ?,23 ... } 

5. The regression multiple correlation coefficients for the n 
equations are computed as follows : 

a,bcde 



ath diag. element 



or -- 



R 



\ Vioi ---' 



y,l23 ... 



Etc. 



Pflcrp ^ 
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6. The standard error of estimate is: 



( 



1 - R< 



7. Universe multiple correlation coefficient 



R = 



/ 



8. 



1 -(i-rO 

Universe standard error of estimate — 



N-l 
N-n 



S = (T (1-R 2 ) 



(— ) 

\ N-l / 



C) The program that has been written assumes the dependent variable 
cells to occupy the last column and/or row. However, a routine 
could be written for any other designation, or since the program 
does not alter the inverse, rows and columns could be interchanged, 
to put the correct elements in the proper locations for this pro- 
gram. 

CODING INFORMATION: 



A) 


Storage 


3 tracks 10 sectors plus externs 


B) 


Linkage 


R (Lo + 309) 
U (Lo ) 


c) 


Input 


None 


D) 


Output 


12.6 - format (See sample proble 




Order: 


h> h l, Rly,2 ... S bl C.R. 

$2, b 2 , R 2y,l ••• S b2 C ' R « 
etc. 

Bn-1, bn-1, Rn-ly,l ... C.R. 

Tab b Tab S bQ C.R. 

C.R. 

A /v. 

R, , S , R , S 


E) 


Constants 


Lo +226 1 @ 29 

Lo + 2ii0 1 @ 2 

Lo + li;8 Floating Point "1" 



F) Timing: Approximately 10 n seconds 

G) Stops: Lo + 308 B.P. k stop at completion of program 
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0000 



00U3 



Initialize 
Lf Matrix 
Lf First Row 
Lo Matrix 



and Lf 6~ ' s 
reg. coeff. 
and Lf Means 



Lo 
Lo 
Lo 
Lo Floating Point 



OOUli 



!_ 



oo55 



Calc. and Float 
N-l 
N-n 



00^6 



i 



onU 



Calc. and Store 

(1-R 2 ) 
^2 



R' 
_S_ 



oii5 



JL 



Calc. 



and Print 

Pi 
b xi 

£xiy,12 . . 
b bxi 



0200 



i + 1 



x. 
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Calculation of Correlation 
and Regression Statistics 



No 



0205 



Yes 



Calc. and print 
^0 y S b0 




0309 i r 



0308 



BP I4 
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JOB NO. 


PROGRAM NO. 

F2-129 


PROGRAM PREPARED BY: [PROGRAM CHECKEO BY: 

Burggrabe j POOL Review 


DATE A 

1-15-60 V 


1 


. ROBLEM: 

External Storage 


TRACK 

62 




PROGRAM INPUT CODES 


0. 



h 
w 


LOCATION 


INSTRUCTION 


a. 


h 
l/l 


CONTENTS 
OF ADDRESS 


NOTES 




OPERATION | ADDRESS 




1 


; 


















/ 


X 















6l 2l3 i 2 
















1 1 3 ,3 


, iXiZj, i i] 




n @ 29 Orde 


r of Inverse * 








l 1 3 ,4 


, ,x , z,[ , , ,] 




N @ 2k No c 


f Records * 








1 |3 i 5 






X 




4 


■ 






1 I 3 1 6 








I 


1 






1 I 3 1 7 


, ,x,z [, , ,] 




Lo Floating 


Point * 








1 l 3 l 8 


■ .I.Z|[. i .] 




Lo Means 


?*~ 








1 I 3 1 9 


, il, z,[, , ,] 




/\| Lo Standai 


■d Deviations * 








1 I 4 1 ° 


i |X , z [ , , ,] 




Lo Inverse 


( D rij Matrix) '"' 








l I 4 1 1 


1 1 1 1 1 1 1 














1 I 4 | 2 


1 1 1 1 1 1 1 














i l 4 i 3 


, ,x, z^, , ,] 




X 


Lo Regress: 


.on Coefficients * | 


| 






1 |4 | 4 
















1 |4 i 5 
















i ,4 i 6 






Note: 


These locations must 








1 |4 i 7 






7, 


be filled ] 


>rior to entering 








1 i 4 1 & 






into subroi 


itine. 4 


I 






l l 4 l 9 








I 


¥ 






l |5 , 
















1 |5 i 1 






X 












1 |5 i 2 
















l |5 i 3 
















1 i 5 | 4 












- - 








1 |5 | 5 




X 












1 |5 i 6 
















1 |5 i 7 
















1 |5 i 8 


[,,,,,,,] 




Temp 










1 |5 i 9 


[,,,;,,,] 




X 


Temp 










1 |6 i 


[..,,.. J 


] 


Temp 


4 


I 






l |6 | 1 


[, , , 1 , , ,] 


i + 1 @ 29 


* 




j 1 1 1 1 1 1 1 _.._ 




1 L 6 1 2 


1 1 1 1 1 1 1 






\ , , , \ , , , 




i .6 , 3 


1 1 1 1 1 1 1 


X 








FORM LP-12 
















CARRIAGE RETU 


RN 


WI..AI M Ul 1 j i :i^4 t X 
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JOB NO. 


PROGRAM NO. 

F2-129 


PROGRAM PREPARED BY: PROGRAM CHECKED BY: 

Burggrabe jPOOL Review 


DATE 

1-15-60 


PROBLEM: A A 

p x1 b xi , Partial Corr. Coeff, S bxi , b Q , & est, R, (Test, R, (f est 


TRACK 

00 




PROGRAM INPUT CODES 


0. 



t- 
m 


LOCATION 


INSTRUCTION 


CL 



h 
i/> 


CONTENTS 
OF ADDRESS 


NOTES 


OPERATION ADDRESS 


1 
1 1 1 1 1 1 1 


1 














1 
1 1 1 I 1 I 1 


1 


X 




_J J 1 1 1 1 1 




, 


,0 ,0 


i i X|Bi 6, 2,3 i3 




n® 29 










[0 i i 


i , X,N. 6 2i3 i3 




n@ 29 










,0 , 2 


i i i M - Oi 2 1 U i0 




1@ 2 




1 I 1 | , 1 1 






l° 1 3 


i i i S | Oi 2 i 2 16 




X 


1 ® 29 




1 






,0 | 4 


i , X|A - 6, 2 |It ,0 






Lo a Matrix 


Inverse 








,0 , 5 


i i i Y . 0, , 6 ,0 




Add [1/V n ] 










,0 | 6 


i i X.Bi 6 2,3 i3 


nft 29 










,0 | 7 


i i i S 1 Oi 2 i 2 16 




X 


1 ft 29 










1° 1 8 


, , X, H 6, 2 , 5 ,U 




Temp 


ctr 








|0 | 9 


i i XiAi 6 2iU i0 




Lo a Matrix 


Inverse 








,1 ,0 


i i , Y . Q, 1,1 ,7 




-Pi/V 




> ,,,■],, , 






,1,1 


i i i Y | Q 1 , 3 ,1 




X 


ii 










|1 I 2 


i , X, B 6, 2, U ,0 




Lo n Matrix 


Inverse 








ll ,3 


i i iYi Q li 2 ,7 




Add l/V,- 










il 1 4 


i i XiBi 6 2,3 i9 




Lo (T ' s 










ll | 5 


i i iY| Q 1|2 i3 




X 


Add \<T x, 1 




i :::;:: : 






,1,6 


i i i Y i °\ 1 i ll i 




J! 




| 








l' l 7 


i i , Y. Q l,Uil 




It 










,1,8 


, , x,A. 6, 2,5,U 




n-1 ® 29 




! 








,1,9 


i i i Y, q 1,1,2 




X 


Add m ] 




: 








,2,0 


i i , Y, Q 1,2 ,2 




II 




1 

! 








|2 i 1 


i i i Yi Q 2i5i2 




II 




i 








|2 | 2 


i i XiB) 6 2ilii3 




Lo reg. Coe 


ff. 








,2 | 3 


i i i y. q 2,3,1 




Xl- 


Add [b n ] 










,2 , 4 


■i i i A, q 2,2,6 




1ft 29 










|2 | 5 


■ i i i Y, q l,2,U 




Add [b^] 










|2 | 6 


| Y 1 2,1 














|2 i 7 


- i i X,B, 6 2,3,8 




X 


Lo Means 




»,,,,,, 






|2 , 8 


i i i Yi 2iOi9 














|2 i 9 


i i X, A, 6, 2,£,li 


_. 


— 


n-1 ^29 










,3,0 


i i i Y, q 2, 2,9 


Add Y 










,3 , 1 


, , X, B. 6, 2,3,7 




X 


Lo Floating 


; Point 


FORM LP-10 














X 


CARRIAGE RETL 


RN 


BO*AL M«»l 


, Jl 364 1 X 



LGP-30 CODING SHEET 



PREPARED FOR: 

LGP-30, RPC-H000 Users' Organization - POOL 


PAGE OF 

2 / 7 


JOB NO. 
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PROGRAM PREPARED BY: (PROGRAM CHECKED BY: 
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DATE 

i-i5-6o 


PROBLEM. A * 

P xi , \^ , Partial Corr. Coeff, S^^. , b Q , & est, R, (Test, R, cT est 


TRACK 

00 


PROGRAM INPUT CODES 


0. 



LOCATION 


INSTRUCTION 


a. 


h 


CONTENTS 
OF ADDRESS 


NOTES 


OPERATION ADDRESS 


1 


/ 








/ 








/ 


X 








1 1 1 1 1 1 1 




,0 ,3 ,2 


i i i y| Oi "Ji7 








1 






|3 i 3 


i , , T, Q Oi 5i8 


; 
/ 












,3 | 4 


, , ,y, qi.i,* 












,3 | 5 


, , ,Y 0,11 6 


X 












,3 | 6 


, , ■*, Q 2,0,7 












I 3 i 7 


_j . .*| ^ 2 i°. 8 . 


/ 
/ 

/ 

/ 
/ 

/ 
i 
i 

i 
i 
i 
i 
i 


— ■ 












,3 | 8 


i i i Y| Q 2,2,7 












,3 , 9 


, , , X, Q 2,2,8 


X — 










I 4 1 ° 


, i i Y, Q 2, Uil 












|4 , t 


i i i Y| 2i ki 2 














,4 , 2 


„ i i i Y| 2, 5i8 












,4 , 3 


i __i i Y| Q 2i^i5> 


X 












l 4 1 4 


, , ,S, 0,2,2,6 




n@ 29 










|4 | 5 


1@ 29 










,4 | 6 


i i X,R, 0, 9,2 £ 


float 


25. OR 








,4 , 7 


i i x i u i °i 9 1 ,0 


X 












|4 i 8 


i i I, Z . 0, 0,2 ,9 














,4 | 9 


i , X, C| 6, 2,5 ,7 


n-1 


' 








,5 , 


i i X, B 6, 2 , 3 ,h 


n @29 










,5 | 1 


. .I,S, 6,2,3,3 


X 


n® 29 










,5 | 2 


i i X,R| Q 9i2 6 


float 


25. OR ' 










l |5 i 3 


, , X, U | 0, 9 , ,0 


i 
i 
i 
i 


— 














|5 i 4 


i i X,Z, Q 0,2 ,9 












|5 | 5 


, , X, C. 6, 2,5 ,8 


X 


n-P 










|5 i 6 


, ,1,0, ^ 2,5 ,3 


Temp -* 










|5 | 7 


i i i R | [ i i 1 3 


I 


2U.0 










|5 | 8 


, i ,U|[ , , ,1 


i 
i 

i 
i 
i 
i 












|5 i 9 


. . ,b| qi,U,8 


X 


in II 










|6 | 


i , ,D|[ , , ,] 




i/v 










L_1 6 „l' 


„„x_„_i_XiH|_..i 2l5 i9 
i . i Xi H| A ?\ ^ 18 


n 


- " ■-—-^r-^r. 








,6 | 2 








,6 | 3 


i i Xi Hi 6 2i5 18 


X 




Ci-R ? ) 




M T 




FORM LP-12 














X 


CARRIAGE RETL 


RN 
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JOB NO. 


PROGRAM NO. 

F2-129 


PROGRAM PREPARED BY: 

Burggrabe 


PROGRAM CHECKED BY: 

POOL Review 


DATE 

1-15-60 


PROBLEM: 

B„, , b r1 , Partial Corr. Coeff , S bxi , b Q , & 


3St 


A A 

, R, (Test, R, <T est 


TRACK 

01 


PROGRAM INPUT CODES 


0. 



I- 


LOCATION 


INSTRUCTION 



h 


CONTENTS 
OF ADDRESS 


NOTES 


OPERATION | ADDRESS 


I 


; 
















/ 


X 








_o, 


1,0 i 


i , X,U j 0, 0,0 ,0 




Acc. -*- M 










,0 | 1 


, I X , M 6 I 2 I^I 7 




(N-l) -* 


(1-R ) 








,0 , 2 


i | X,Y 0, 0,0 ,0 




change signs -*• 


-1 + R^ 








,0 | 3 


i | , A . 0, 1 1 k ,8 




X 


tin n ^. 


A O 

R^ 








,0 | 4 


, , XiH | 6, 2,6 ,0 














,0 | 5 


i i X,B| 6, 2,5 i7 




N-l 










,0 ,6 


i i i A | o, l , h |8 


T 


1 










1° 1 1 


, | X,H . 6, 2,5 $ 


X 


-> 


n 








,0 | 8 


i i XiU | Oi OiO i0 




A O 










|0 | 9 


i i X i M i 6 i 2 i^ i 8 




l-IT 


(1-R) C N 








,1,0 


, , X,R - 0, 0,0 ,0 




V~ 


N-l 








,1,1 


, , X,U, 0, 0,0 ,0 




X 












1 1 i 2 


, , ,M,r , , ,1 




<Tj 










|1 l 3 


i , X,H . 6, 2,5 ,8 




A 

s 










i1 I * 


i i X|E - 0, 0,0 ,0 




Exit 


F.P. 








|1 i 5 


, , iR|[ , i .] 




X 


2U.0 










,1,6 


, , ,u,[ , , ,] 














,1 ,7 


, , ,p,[ , I ,1 




-P/v 










,1,8 


i i X.Mj 6 2,5 .9 




n 










,1,9 


, | x,y q o,o ,0 




X 


Change Sigr 


s 


/ 






,2 ,0 


. , X,P 0, 0,0 ,0 




Print B t 










|2 i 1 


( , X,U 0, 0,0 ,0 




Acc — >- M 










|2 i 2 


i i t m | r i i 




<ry 










|2 | 3 


, , ,d][ , , ,] 




^L 


^xi 










,2 , 4 


, , ,H,[ , , ,] 




b . 










(2 i 5 


i , X|P. 0, 0,0 ,0 




XI 










|2 i 6 


i , , B , 0, 1 1 h ,8 




in it 










|2 i 7 


, , ,D,[ , , ,] 


X 


i/v_. 










|2 | 8 


, , X,H. 6, 2,5 ,7 




(i-Ri I 










|2 i 9 


i i XiM| 6 2,5 i7 


Pi -Vi 




I I 1 I I I I 






,3 , 


, , X,U, Q 0,0 ,0 








1 1 1 I 1 1 1 






|3 i 1 


i , i m | r i i ii 




X 


-PiAi 


-Piv'Pvi 


FORM LP-IO 














X 


CARRIAGE RETt. 


JRN 


nor.L »«.!(. 4UI4K 
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JOB NO. 


PROGRAM NO. 

F2-129 


PROGRAM PREPARED BY: 

Burggrabe 


PROGRAM CHECKED BY: 

POOL Review 


DATE 

1-15-60 






PROBLEM. ^ ^ 

Pxi> b xi> Partial Corr. Coeff, S bxi , b , & est, R, crest, R, ..-Test 


1 RACK 

01 






PROGRAM INPUT CODES 


0. 


h 


LOCATION 


INSTRUCTION 


a. 


h 

l/l 


CONTENTS 
OF ADDRESS 


NOTES 




OPERATION | AODRESS 




1 
1 I 1 I 1 1 1 


1 












--- - - 




1 1 1 | 1 1 1 


1 


X 















0,1 ,3 ,2 


i i I,T| 0, 0,0 |0 












I 1 3 ,3 


, , X,H ( 6, 2,5,6 


— 


2 

rij .xxx 


Partial 








i 1 3 ,4 


i i XiR, 0, 0,0 ,0 




V" 


rij-xx? 


c 








1 |3 i 5 


, , X, P 0, , ,0 




X 












I I 3 I 6 


. , , B , °^^P 






IITJt 










I I 3 i 7 


i i Xi s I 6, 2 1 5 16 


r-2 










I i 3 i 8 


i i X|D| 6, 2,5 ,7 


(1-R 2) 




— 








1 i 3 1 9 


i i x, D ! 6, 2,5 £ 




£A w 










1 I 4 1 ° 


, , iD.[ , , ,] 




6 














i J 4 1 ' 


D [ ] 

1 1 1 r 1 1 l J 


XI 


fr-^/Ntr (1-R-) 








_L L 4 1 2 


i i XiR 1 0i 0,0 ,0 


T 


A. X 


VcT. 








i i 4 i 3 


i i X,U | 0, 0,0 ,0 


X 


Ace -*■ M 


i'j 


- 








1 I 4 1 4 


x. i_XiM^ 6 L 2i5 i8 




A 

<f~ est 










l I 4 1 5 


1 1 XiP| Ol O'O :0 




Print 


(s,.) 








1 L 4 1 6 


i i Xi E , 0i 0,0 ,0 






bxi 








1 l 4 i 7 


i i X,P, 1, 6,0 ,0 


I 


y 








, ,o p 9 o,o ,o i 


1 


l i 4 l » 


V 9 0,0 0, 0,0 ,2 




F.P. "1" &| 


delay 








1 l 4 l 9 


, i ,Bj 0,1,1 ,7 




Add [-PiA n ] 










l |5 , 


i i X, A | 6: 2,3 ,3 




n@ 29 










1 |5 i t 


, , J 0,1,1 ,7 




X 












l |5 i 2 


i i i Y | 0, 1 1 3 1 














l |5 i 3 


, , , B 0, 1 , 2 ,3 




^I^xlL 










1 |5 i 4 


i i i A | 0i 2 , 2 ,6 




la 2? 




- - - - 








1 |5 | 5 


■ , ,1,0,1,23 


7 


X 












l |5 i 6 


, , ,T, 0, 1,U P 












1 |5 i 7 


i i i Y | 0, 1 , h 1 














i 1 5 ,8 


i , i B | 0, 1 , 2 <1+ 




, H[b .] 










i |5 i 9 


i i , A | 0: 2 , 2 ,6 


' 


7 


xr 
1® 29 










1 ib , 


i , iY . 0. 1,2 h 


----- 










- 1 i L 1 J_ 1 1 




I ,6,1 


_^ , B . O 12 7 


D[lAj ] 
(n+1) -5 29 




1 J 1 III i 




I ,6,2 


, , X.A . 6. 2 .6 1 




1 - 
1 1 1 | 1 1 ! 




i ,6,3 


, , ,Y , 0, 1 ,2 ' 7 


i7\ 


f 


1 
'OHM LP-t2 






, , 

/ 

/ \ 
t— —J 


CARRIAGE RETU 


RN 




X 
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PROBLEM: A A 

B xi> b xi> Partial Corr. Coeff , S bxi , b , & est, R, <T est, R, <T est 


02 


PROGRAM INPUT CODES 


a. 


u 


LOCATION 


INSTRUCTION ^ 


CONTENTS 
OF ADDRESS 


NOTES 


OPERATION ADDRESS J 


1 


; 
















X 










,0 ,0 


i 1I1B | 6.2.5 3 ' 


Ctr 










,0 i 1 


i , ,A 0,2 ,2 6 ' 


1 @ 29 










|0 | 2 


i ,X,H 6,2,5 3 ' 


Ctr 










,0 | 3 


i ,x,s , 6,2,5 h ' 


:x 


n - 1 @ 29 










,0 |4 


i , iT . 0,1,1 $ ' 












,0 , 5 


, ,X,C ,6,2,5 3 / 


Ctr 










,0 ,6 


, iIiC j 6,2,5 6 ' 


Z Area 










,0 ,7 


, . .R,[ , . ,1 ' 


:z 


2U.0 










,0 ,8 


i . iU.[ , , ,] ' 












,0 ,9 


, , ,p,[ , . ,1 ' 


X. 

i 










,1,0 


, . ,M ,[ , , ,] ' 


^, h ^ 




1 1 1 I 1 1 1 






,1,1 


, , X,A . 6, 2,5 ,6 


„x 


Z 










|1 |2 


. , X,H 6, 2,5 ,6 


z 


1 








1 1 | 3 


X E 
I . i | i . . 


Exit F.P. 










il l 4 


i i , B 0, 2 , ,9 


' P[ x\ ] 










,1 i 5 


i A , °i 2 l 2 , 6 


'X 


1@ 29 










,1,6 


, , , Y 0, 2 , ,9 












,1 ,7 


i . iB | 0, 2,1 ,0 












,1,8 


i i , A , 0, 2 , 2 ,6 


' 1@ 29 










,1,9 


i i , 1 . 0, 2 , 1 ,0 


'X 












,2 ,0 


. . X,B 6, 2,5 ,3 


' ctr 










i |2 i 1 


, , , A , 0, 2 , 2 ,6 


' 1@ 29 










|2 |2 


i i XiH. 6 2.5 ,3 


' ctr 










,2 | 3 


. . x.s. 6 2,5 iU 


-X 


(n-l)@ 29 










,2 , 4 


, . , T . q 2,0 ,7 












,2 i 5 


, , X, P. 2, U,0,0 


' tab 










|2 , 6 


, , X, Z , 0, 0, ,1 


' 1@ 29 & c 


el ay 








|2 i 7 


, , ,R,r , , ,] 


lX 


2U.0 










|2 i 8 


, , ,u,r , , ,] 












I |2 |9 


, , ,B,r , , ,] 


' Y 










I |3 i 


i i X, S| 6 2,5 16 


' Z 










i |3 i 1 


, , ,H,r. , ,i 


'X 


b o 




FORM UP- 10 












X 


CARRIAGE RETl 


JRN 


KOYAl M«BC£. j)}MI X 



rnNniTinNAi stop r.onF 
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PROGRAM NO. 

F2-129 


PROGRAM PREPAREO BY 

Burggrabe 


PROGRAM CHECKED BY: 

POOL Review 


DATE 

1-15-60 


PROBLEM: , » 

ft, n - , \i , Partial Corr. Coeff, Sb xi , b n , & est, R, (T est, R, rf est 


TRACK 

02 


PROGRAM INPUT CODES 


u. 


i- 


LOCATION 


INSTRUCTION 


a. 


1- 

V) 


CONTENTS 
OF ADDRESS 


NOTES 


OPERATION 1 ADDRESS 


I 


/ 
















/ 


X 












0,2 , 3 , 2 


i i X, P. Q OiOiO 










1 1 1 1 1 1 1 




1 |3 i 3 


i , x,B 6, 2,5 $ 


— 


N 








l |3 | 4 


, , X,R 0,0,0,0 


-r -*■ 


•VF 






l 1 3 | 5 


i i IiH, 6, 2,5 $ 




y. 










1 | 3 | 6 


i , X,B 6, 2,5 ,8 




A 

S est 








I I 3 i 7 


i i X|D - 6, 2,5 $ 




-/¥ 








1 i 3 |8 


, i X, E | 0, , ,0 




Exit F.P. 








1 i 3 i 9 


, , X, P .2, U|0 ,0 




Xltab 




,000 ( ( o r 1 


1 


1 I 4 |0 


2000 0000 
I I i 1 i i i 




1® 2 and d 


slay 






1 I 4 I ' 


, , ,R|[ , , , ] 












1 I 4 1 2 


, , ,u,r ,i,i 












i l 4 i 3 


, , X,P 0, 0,0 ,0 




X 


Sb 








1 |4 i 4 


i i X, B j 6, 2,5 ,9 




1-R2 








l I 4 i 5 


i , i S | Q 1 , U ,8 




IITJI 








i i 4 i 6 


, , X,Y 0, 0,0 ,0 




Change sign 


s R^ 






1 I 4 i 7 


, X, R 0, , ,0 




X 


v~ 


R 






l |4 i 8 


i , X.Hj 6, 2,5 ,7 




R 


1 






1 (4 i 9 


, , X , B 6 , 2 ,£ ? 




1-R 2 


1 






l |5 , 


i , X, R | 0, , ,0 




V~ 








1 |5 i 1 


, , X, U . 0, , ,0 




X 


Ace -*■ M 








1 |5 i 2 


, , ,M,[ , , , ] 




<TJ — 


S est 






l |5 i 3 


i i X , H | 6 , 2 i^ ? 




S est 








1 |5 l 4 


i , X, E | 0, , ,0 












1 |5 t 5 


, i XiP, l, 6,0 ,0 




X 


cr. 








1 |5 i 6 


i , X,Z | 0, 0,0 ,0 












1 |5 i 7 


, , X,P, 1, 6,0 ,0 




cr. 








l |5 i 8 


, i ,R|[ , , , 1 












1 |5 i 9 


, , ,u ,r , , , i 




X 










1 |6 i 


i i X,B, 6, 2,5 ,7 


__ 


R 








I |6 i 1 


| I | P 1 0,0^^ 


Print 








1 16 t 2 


, , X, B , 6, 2 , 5 ,9 




S est 








i 1 6 ,3 


1 
i i X, P, 0, 0,0 ,0 


X 


Print 




FORM LP- 12 














X 


CARRIAGE RETU 


RN 
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• 


OB NO. 


PROGRAM NO. PROGRAM PREPARED BY: 

F2-129 Burggrabe 


PROGRAM CHECKEO BY: 

POOL Review 


1-15-60 


F 


' R ° BLEM: 3^, *xi> P«tial Corr ' Coeff > Sbxi ' bo ' & eSt ' R ' ****' R ' ^ 6St 


03 




PROGRAM INPUT CODES 


CL 



h 


LOCATION 


INSTRUCTION q 


CONTENTS 
OF ADDRESS 


NOTES 




OPERATION ADDRESS Jj 




1 




1 


















X 










Oj. 


3,0)0 


, , X,B, 6, 2,6,0 ' 


R^ 












|0 |1 


, x,R 0, 0,0,0 l 














,0 ,2 


i , x, p q o, o ,o ' 














,0 ,3 


, , X,B 6 2,5,(i ' 


.X. 


5 est 




• 








,0 | 4 


, , x,p q 0,0,0 ' 








1 | | | | | | 






,0 ,5 


s , x,E, q 0,0,0 / 














i° I 6 


| , Xi Pi 1 6,0,0 ' 


cr. 












,0 ,7 


i , x, z. q o, o,o ' 


.X. 


delay 












|0 |8 


i i x, z . q U, o ,o ' 


Stop 


BP h 










,0 ,9 


, , . u,r . , .1 ' 


Exit 












,1 ,0 








• 








|1 i' 




.X 












|1 |2 
















il |3 
















,1 i 4 


1 1 ! 1 > I 1 














,1,5 




'.X 






A 








i il |6 








9 








i .1 i7 
















i 1 1 i 8 
















i ,1 |9 




lX 














i |2 ,0 
















i |2 i 1 
















i i2 ,2 
















1 |2 | 3 




lX 














i i2 ,4 
















l |2 i 5 
















i |2 ,6 
















l |2 |7 




'X 






• 








1 |2 i 8 


1 1 1 1 1 l 1 






1 1 1 ^_J 1 L 






I |2 |9 


1 














i .3 ,0 
















i • 3 i 1 




'X 




■ 




FORM LP- 10 












X 


CARRIAGE RET 


URN 





CONDITIONAL STOP CODE 
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TITLE: Calculation of Y 1 and (y-j 1 ) given a set of coefficients 
(Bo, Bi, ... B n ) and a set of records (l, X]_, X 2 ... X n ) 

AUTHOR : Allen G. Renz 

Compumatix, Incorporated 

DATE : October 8, 1959 

PURPOSE : To read in a set of records, each record containing values 
for X]_, X 2 , . . . X n and Y, and calculate a predicted value 
of Y, Y 1 , using a set of coefficients stored in the machine. 
Also calculate (j-j 1 ) and print out Y, Y 1 , and (j-j 1 ) . 

RESTRICTIONS : 

A) Normal restrictions of 11.2 and 12.1 

B) Coefficients must be stored in machine prior to 
entering subroutine in order, (i.e., b 'b -^ 

b xn^ in fixed Point @ q of variables. ' ' 

C) The number of variables is limited only by machine 
storage available. 

CODING INFORMATION : 

A) Storage 

1) Program - Sh sectors 

2) External storage - 6232 to 6263 (See attached sheet) 

3) 11.2 data input - 0300-0563 
h) 12.1 data output - O6OO-O763 

B) No calling sequence is required. 
Linkage: R (Lo + 50) 

U (Lo + 00) 

C) Input is in 11.2 data input format in the following 
sequence: 
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7 



First record 

P ± qq (Lo record) ' 1' X x ' X 2 ' . . . X n l Y» -0000000" 
Following Records: 

P ± qq (Lo record +1) 1 X x « X 2 ' ... X n I Y 1 -0000000" 

D) Output: Output is in 12.1 fixed point data output 
format. Three colums are printed out: 

T x Y 1 and (y-y 1 ) at ZQ of variables (X 1 ... X n ) 

E) Location of constants: 
Lo + 36 18 29 

Lo + 51 XZ 0300 
Lo + 52 XZ 6363 
Lo + 53 Temp storage 

F) Timing: 5 sec/l5 X's/ record + input and print time 

G) Program stops: 

Lo + ii5 B.P. k after print out 
Lo + l\9 B.P. 8 after completion 

H) Storage Requirement External) 

6232 - Two times the data q @ q = 29 

6233 - n, no. of variables @ q = 29 
623U - N, no. of records @ q = 29 
6235 - Lo of the record 

62U3 - Lo of the coefficients 
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7 



Calculation of Y 1 and (y-y 1 ) 



Initialize 



Data Input 



I EXIT ] 



Calculate 
2 x b = y» 





Test for 
End of Data 
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JOB NO. 


PROGRAM NO. 

F2-129 


PROGRAM PREPARED BY: 

Burggrabe 


PROGRAM CHECKED BY: 

POOL Review 


DATE 

10-6-59 


PROBLEM: 

Calc. & Print Y, Y' (Y-Y 1 ) Fixed Point 


TRACK 

00 


PROGRAM INPUT COOES 


&. 


1- 

1/1 


LOCATION 


INSTRUCTION 


Cl 



CONTENTS 
OF ADDRESS 


NOTES 


OPERATION ADORESS 




1 


1 


















X 








_3 o,o ,o 


i i X, B| 6 2,3,5 




Lo Record 










,0 | 1 


i i ,y, o on 18 














,0 ,2 


i i i Y, Q 0i3i0 














,0 ,3 


i i XI A. 6, 2,3,3 




X 


n@ 29 


no variables 








,0 | 4 


1 . i Y , 9 0, 3 ,1 






Lf Record 








,0 ,5 


i . . Y i 9 °i^. 3 














,0 ,6 


i i X.B. 6, 2, U,3 




Lo of coeffic 


ients 








,0 ,7 


i i i Y, q 0,1,7 




X 












,0 ,« 


i i i B. Q 0|£|1 




XZ0300 










,0 ,9 


, , X, A, 6, 2,3,2 




Q of printc 


tit 








,1,0 


i . i y, q o, U,o 














,1,1 


i , X, C, 6 2,6,0 




X 


Ace -»0 










i» |2 


, i X S ( 6, 2, 3,U 




M@ 29 


Number of records 








1 1 i 3 


i , x c i ^ 2 i 6 i° 




ctr 










i1 1 * 


i i X C. 6. 2,5,8 




2 storage 










1 1 , 5 


i i X R i 9 3, 0,8 




X 












,1 ,6 


i , X Uj Q 3,0,0 














i' |7 


, , ,B,[, , , ] 




b . 










,1,8 


. i .M,[, , , ] 




XI 
X. 










,1,9 


i i X A 6. 2, 5,8 




X 


2(b-X), 










,2 ,0 


, , X C| 6 2i 5,8 




Z(b-X)i + 1 










i |2 i 1 


i i , B i 9 0|1.7 




B M b xi 










|2 |2 


i , , A, q o, 3,6 




1 @ 29 










,2 | 3 


, . . Y, 0, 1,7 




X 












,2 , 4 


i i , Bi 0, 1,8 




M [ ] x ± 










|2 i 5 


i , , A q 0, 3 ,6 




1® 29 










,2 | 6 


i , iT 1 Q 0,1 ,8 














|2 i 7 


, i iEi 0,5,2 




X 


xz 6363 










|2 i 8 


i , , S | Q 0,5,3 




T, r record 










,2 | 9 


i i i T 1 Q 0,li7 




i 










,3 | 


, , ,B |L ,1,1 




10 a 










|3 i 1 


, , ,M,[ , , , ] 




X 


y ' q 


y@ 2q 


FORM LP-10 










L -" I 


/ 


X 


CARRIAGE RETL 
CONDITIONAL S 


RN 

TOP CODE 


hoial H'lu, J13U41 * 
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JOB NO. 


PROGRAM NO. 

F2-129 


PROGRAM PREPARED BY: 

Burggrabe 


PROGRAM CHECKED BY: 

POOL Review 


DATE 

10-6-59 i 


PROBLEM: 

Calc. & Print Y, Y» (Y-Y 1 ) Fixed Point 


TRACK 

00 


PROGRAM INPUT CODES 


0. 


h 
VI 


LOCATION 


INSTRUCTION 




1- 

10 


CONTENTS 
OF ADDRESS 


NOTES 


OPERATION 1 ADDRES5 


1 


/ 
















/ 


X 








_0_, 


,3 ,2 


i |X,H , 6,2,5 ,7 




Y 


7 








|3 i 3 


i iX,S 6,2,5-3 






Y' 

Y-Y' 


U Q 








,3 | 4 


i iXiC | 6i2i5 9 












|3 i 5 


i 1X1P | I161O o 




X 


cr. 




1 


m 




,3 | 6 


i iX,Z , 0,0,0 I 




delay & 1 


9 29 








,3 | 7 


i |X,B | 6,2,5 ,7 




Lo print a 


rea 








,3 | 8 


i i XiR | 0i6,0 5 


12.1a 










,3 | 9 


1 
i i XiU | Oi 6i0 




Xl 












,4 | 


i 1X1Z |[ i i i 1 




NNQQ 




1 1 1 I 1 1 1 






,4 , 1 


i i XiB , 6, 2,3 5 




Lo record 










,4 | 2 


i i iY | OiOil fi 














,4 ., 3 


, ,X,B 6,2,1*3 




X 


Lo coeff. 


_ „ _ _ | 








|4 | 4 


1 
j. I lY i DlQjIjZ 














|4 | 5 


i i X i Z °i ^i° i° 




BP h after 


each record 








|4 i 6 


i i XiB | 6i 2,6 i0 














|4 i 7 


i i , A , 0, 1 3 16 




X 


1 @ 29 










|4 i 8 


i i it | oi on a 






)| 








|4 i 9 


i i x, Z | 0, 8 , ,0 




BP 8 


At completion 








,5 | 


, , ,u,[ i , , 1 




Exit 










|5 i 1 


i i X ,Z 0, 3,0 ,0 




X 


for print 


3ode (12.1a) 








,5 | 2 


i i X, z |6 3 6 \ 




mask 










|5 i 3 


, , X,Z | 0, 0,0 ,0 




Lf Record 










|5 | 4 


1 1 1 1 1 1 1 














i5 i 5 


1 1 1 1 1 1 1 




X 












i5 i 6 


1 1 1 1 1 1 1 














|5 i 7 


1 1 1 1 1 1 1 














|5 i 8 


1 1 1 1 1 1 1 














|5 i 9 


1 1 1 1 1 1 1 




X 












,6 i 


1 1 1 1 1 1 1 


— 














,6 | \ 


1 1 1 1 1 1 I 














,6 | 2 


1 1 1 1 1 1 1 


— 














,6 | 3 


I 1 1 1 1 ! 1 


X 






FORM LP-12 














X 


CARRIAGE RETL 


RN 


"'■'•I «"• ' J»5!>4» X 



CONDITIONAL STOP CODE 



